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ABSTRACT 

We  present  an  efficient  algorithm  for  planning  the  motion 
of  a  convex  polygonal  body  B  in  two-dimensional  space  bounded 
by  a  collection  of  polygonal  obstacles.  Our  algorithm  extends 
and  combines  the  techniques  of  Leven  and  Sharir  and  of  Sifrony 
and  Sharir  used  for  the  case  in  which  B  is  a  line  segment  (a 
"ladder").  It  also  makes  use  of  the  results  of  Kedem  and  Sharir 
on  the  planning  of  translational  motion  of  B  amidst  polygonal 
obstacles,  and  of  a  recent  result  of  Leven  and  Sharir  on  the 
number  of  free  critical  contacts  of  B  with  such  polygonal  obsta- 
cles. The  algorithm  runs  in  time  O  {kn\j{kn)  log  kn),  where  k  is 
the  number  of  sides  of  B,  n  is  the  number  of  obstacle  edges,  and 
\j(<7)  is  an  almost  linear  function  of  q  yielding  the  maximal 
number  of  connected  portions  of  q  continuous  functions  which 
compose  the  graph  of  their  lower  envelope,  where  it  is  assumed 
that  each  pair  of  these  functions  intersect  in  at  most  s  points. 


1.  Introduction 

Let  B  he  a  convex  polygonal  object  having  k  vertices  and  edges,  free  to 
move  (translate  and  rotate)  in  an  open  two-dimensional  space  V  bounded  by  a 
collection  of  polygonal  obstacles  ("walls")  having  altogether  n  comers.  The 
problem  studied  in  this  paper  is  to  plan  automatically  a  continuous  obstacle- 
avoiding  motion  of  5  between  any  two  specified  initial  and  final  placements; 
cf.  Fig.  l.L 
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This  problem  has  been  firs:  considered  by  Schwartz  and  Sharir  [SS],  who 
present  an  0{n^)  algorithin  for  its  solution  (which  applies  to  non-convex 
polygonal  moving  objects  as  well).  Since  then  this  classical  "piano-movers" 
problem  had  been  studied  exier^ively,  and  several  efficient  algorithms  had 
since  been  developed  for  certain  special  cases  of  it  [BZ],  [LSI],  [LS2],  [KS], 
[KLPS],  [OY],  [OSYl],  [0SY2],  [SiS]. 

The  algorithm  developed  in  this  paper  is  based  on  a  generalization  and 
combination  of  the  algorithms  developed  by  Leven  and  Sharir  [LSI]  and  by 
Sifrony  and  Sharir  [SiS]  for  the  case  in  which  B  is  a  line  segment.  The 
Leven-Sharir  algorithm  partitions  the  (3-dimensional)  space  FP  of  free  place- 
ments of  B  (also  known  as  the  free  configuration  space  of  B)  into  simple, 
openly  disjoint  and  connected  cells,  and  then  determines  the  adjacency 
between  these  cells.  This  yields  aii  abstract  representation  of  FP  by  a  connec- 
tivity graph  whose  nodes  are  ttiese  cells  and  whose  edges  connect  pairs  of 
adjacent  cells.  Once  the  cells  containing  the  specified  initial  and  final  place- 
ments of  B  are  determined,  the  motion-planning  problem  is  then  reduced  to  a 
simple  graph  searching. 

The  second  algorithm  [SiS]  also  reduces  the  problem  to  a  combinatorial 
graph  searching,  but  uses  a  different  graph,  called  the  vertex  graph,  whose 
nodes  are  the  comers  of  FP,  and  whose  edges  connect  pairs  of  comers  that 
are  adjacent  along  edges  of  FP,  or  along  some  additional  auxiliary  arcs  in  FP. 

Our  algorithm  constructs  an  "intermediate"  kind  of  graph,  which  we  call 
an  edge  graph;  its  nodes  are  edges  of  the  boundary  of  FP  and  its  edges  con- 
nect pairs  of  adjacent  FP-edges  (in  a  sense  to  be  defined  more  precisely 
below).  Our  algorithm  begins  (as  in  [LSI])  by  restricting  the  motion  of  B  to 
be  purely  translational  at  some  fixed  orientation  9 .  This  motion  has  only  two 
degrees  of  freedom,  so  that  it  is  easier  to  calculate  its  associated  restricted  7- 
D  space  FPq  of  free  placements  (a  task  which  has  already  been  carried  out  in 
[KS];  cf.  also  [BZ],  [KLPS],  [LS2]),  and  to  represent  it  as  the  union  of  2-D 
polygonal  regions  having  simple  shape.  Each  such  region  can  be  given  a 
discrete  combinatorial  labeling  that  does  not  depend  continuously  on  6. 
Roughly  speaking,  we  represent  FP^  by  a  graph  VG^,  whose  nodes  are  the 
comers  of  FPq,  and  whose  edges  connect  pairs  of  adjacent  comers.  The 
nodes  of  VG^  are  given  discrete  combinatorial  labels  (that  do  not  depend 
continuously  on  6). 

Next  we  observe  that  this  combinatorial  description  of  FPq  will  not 
change  as  9  varies  slightly,  unless  9  is  one  of  finitely  many  critical  orienta- 
tions, at  which  some  critical  condition,  which  affects  the  combinatorial  stmc- 
ture  of  FPe ,  occurs. 
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As  it  turns  out,  the  most  complex  of  these  criticaJ  orientations  are  those 
at  which  the  object  B  makes  three  simultaneous  contacts  with  the  obstacles, 
without  penetrating  into  any  obstacle.  If  fl  is  a  line  segment  (a  "ladder"), 
then  it  is  shown  in  [LSI]  that  the  total  munber  of  such  critical  placements  of 
B  is  0{n^),  which  consequently  leads  to  an  O(n^logn)  algorithm  for  the 
desired  motion  planning.  If  fi  is  a  convex  it-gon  which  is  free  only  to  translate 
in  V  but  not  to  rotate,  then  the  motion  planning  problem  becomes  simpler 
and  can  be  accomplished  in  time  0{kn  log  kn)  [LS2],  [KS],  [KLPS],  [BZ]. 
This  follows  from  the  property,  which  is  proved  in  [KS]  and  will  be  used 
below,  that  the  number  of  placements  of  B  (all  having  the  same  given  orien- 
tation) at  which  it  simultaneously  touches  two  obstacles,  without  penetrating 
mto  any  obstacle,  is  0{n)  (provided  that  B  and  the  obstacles  are  in  "general 
position";  cf.  [KS]  and  below).  If  B  is  also  allowed  to  rotate,  then  the 
corresponding  critical  orientations  are  much  harder  to  analyze.  Since  each 
contact  of  B  with  the  walls  is  a  contact  of  either  a  corner  of  B  with  a  wall 
edge  or  an  edge  of  B  with  a  wall  comer,  a  crude  and  straightforward  upper 
bound  on  the  number  of  these  critical  placements  of  triple  contact  of  B  is 
0{{kny).  Moreover,  if  B  is  nonconvex,  then  there  are  cases  where  the 
number  of  these  critical  placements  of  B  is  indeed  Cl{{kny).  However,  a 
recent  result  of  Leven  and  Sharir  [LS3]  shows  that  if  B  is  convex,  then  the 
number  of  these  critical  orientations  is  only  0{kn\j{kn)),  where  \j{q)  is  the 
maximal  number  of  connected  portions  of  the  graphs  of  q  continuous  func- 
tions which  compose  the  graph  of  their  lower  envelope,  where  it  is  assumed 
that  each  pair  of  these  functions  intersect  in  at  most  s  points.  It  is  shown  in 
[Sz]  that  Xj(n)  =  0{n  \og*n)   (where  \og*n  is  the  length  of  the  smallest 

exponential  tower  2^  exceeding  or  equal  to  n).    A  better  asymptotic 

bound  is  given  in  [HS]  for  the  case  s  =  3  and  in  [Sh]  for  larger  values  of  s. 
These  better  bounds  are  roughly  of  the  form  0(na(n)'^^'^^"^'  )),  where  a(n) 
is  the  functional  inverse  of  Ackermann's  function,  and  is  extremely  slowly 
growing.  In  short,  for  a  fixed  s,  Xj(rz)  is  nearly  linear  in  n,  although,  as  is 
shown  in  [HS],  [Sh2]  one  has  X,(n)  =  a(na\-^'-'^'>'^i{n)),  so  that  X,(n)  is 
super  linear  in  n  for  j>3. 

Using  these  bounds  together  with  the  techniques  of  [LSI],  [SiS],  [KLPS] 
and  [KS],  we  next  extend  each  node  in  VG^,  by  varying  0,  into  a  node  which 
represents  an  edge  of  FP,  and  then  construct  an  edge  graph  EG  which 
represents  adjacency  of  such  edges  along  the  boundary  of  FP.  All  this  finally 
yields  a  motion  planning  algorithm  for  a  convex  polygonal  object  B  which 
nms  in  time  0{kn\j{kn)  log  kn).  This  algorithm  is  being  implemented  on  an 
IBM  robot  RS/2,  and  its  experimental  results  will  be  described  in  a  forthcom- 
ing paper. 
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The  paper  is  organized  as  follows.  In  section  2  we  describe  the  recursive 
discrete  representation  of  FF  by  its  associated  edge-graph  EG,  as  outlined 
above.  In  section  3  we  discuss  ihe  algorithmic  details  involved  in  the  con- 
struction of  EG,  and  in  its  use  foi  3:ctual  motion  planning. 

2.  Discrete  Recursive  Representation  of  the  Free  Configuration  Space. 

Let  B  be  a  bounded  convex  ^-sided  polygonal  body,  whose  interior  is 
nonempty.  B  is  free  to  move  (translate  and  rotate)  in  a  bounded  two  dimen- 
sional open  region  V  having  a  polygonal  boundary  with  n  corners  altogether. 
We  assume  that  this  boundary  can  be  partitioned  into  a  collection  of  pairwise 
disjoint  simple  closed  polygonal  curves  (which  we  call  "walls");  the  "wall 
region"  V  can  then  be  partitioned  into  a  collection  of  convex  polygonal 
regions  having  pairwise  disjoint  (and  nonempty)  interiors.  Note  that  these 
assumptions  exclude  degenerate  configurations  in  which  a  wall  is  just  an  iso- 
lated comer  or  an  isolated  segment,  or  in  which  a  wall  comer  is  adjacent  to 
more  than  two  wall  edges. 

Let  P  be  an  arbitrary  fixed  reference  point  in  the  interior  of  B,  and  let  Q 
be  an  arbitrary  comer  of  B.  Each  placement  Z  of  fi  in  the  plane  can  be 
represented  by  the  parameters  (AT,  9),  where  X  is  the  position  of  P,  and 
where  0  is  the  orientation  of  the  vector  PQ. 

As  in  [SS],  we  define  a  free  placement  (X,Q)  of  B  to  be  a  placement  at 
which  B  is  fully  contained  within  V;  a  semi  free  placement  {X,  9)  of  B  is 
defined  to  be  a  placement  at  which  B  may  touch  some  walls,  but  not 
penetrate  into  the  interior  of  the  wall  region  V .  The  set  FP  of  all  free  place- 
ments of  5  is  an  open  three  dimensional  manifold,  and  the  set  SFP  of  all  semi 
free  placements  is  closed. 

As  in  [LS3],  we  assume  that  the  moving  object  B  and  the  obstacles  are  in 
general  position.  Roughly  speaking,  this  means  that  (i)  there  does  not  exist  a 
placement  of  B  in  which  it  meets  four  independent  constraints  involving  con- 
tacts with  obstacles;  and  (ii)  there  do  not  exist  two  placements  of  B  with  the 
same  orientation  such  that  B  meets  at  each  of  them  three  independent  con- 
straints involving  contacts  with  obstacles  (cf.  [LS3]  for  more  detail).  We  also 
assume  for  simplicity  that  no  wall  edge  is  horizontal;  this  can  always  be 
enforced  by  an  appropriate  rotation  of  V. 

Since  the  motion  of  B  has  three  degrees  of  freedom,  we  first  analyze,  as 
in  [LSI],  only  purely  translational  motion  of  B  (involving  just  two  degrees  of 
freedom),  and  only  then  treat  the  case  of  general  motion  of  B,  including  rota- 
tion. This  will  enable  us  to  obtain  recursively  a  combinatorial  representation 
of  each  cross-section  of  FP  at  a  fixed  9,  from  which  we  will  then  construct  a 
certain  discrete  graph  which  represents  the  entire  space  FP  in  a  "connectivity 


preserving"  manner,  and  which  allows  us  to  reduce  the  motion  planning  prob- 
lem to  a  discrete  problem  of  path  searching  through  that  graph.  (Our 
method  can  be  regarded  as  a  hybrid  of  the  two  techniques  presented  in  [LSI] 
and  [SiS],  in  a  sense  that  will  become  clear  later  on.  An  attempt  at  direct 
generalization  of  the  technique  of  [LSI]  has  led  to  certain  technical  difficul- 
ties that  we  have  not  been  able  to  overcome) 

2.1.  The  case  of  translatlonal  motion  of  B. 

Definition  2.1:  (a)  ([LS3])  A  (potential)  contact  pair  O  is  a  pair  (W,S)  such 
that  either  W  is  a  wall  edge  and  5  is  a  corner  of  5,  or  W  is  a  wall  comer  and  5 
is  a  side  of  S.  In  the  first  case  we  call  the  pair  a  contact  pair  of  type  I  and  in 
the  second  case  a  contact  pair  of  type  n. 

(b)  A  contact  pair  of  type  III  is  a  pair  O  =  {W,S)  where  W  is  a  wall  comer 
and  5  is  a  comer  of  B. 

(c)  An  actual  obstacle  contact  (i.e.  a  contact  of  B  with  an  obstacle)  is  said  to 
involve  the  contact  pair  O  =  {W,S)  if  this  contact  is  of  a  point  on  5  against  a 
point  of  W,  and  furthermore  if  this  contact  is  locally  free.  i.e.  the  inner  angle 
of  5  at  5  lies  entirely  on  the  exterior  side  of  W  if  5  is  a  comer  of  B,  and  the 
entire  angle  within  the  wall  region  V  at  W  lies  exterior  to  5  if  W  is  a  wall 
comer. 

(d)  (cf .  [KS])  Let  A  be  one  of  the  convex  polygonal  obstacles  into  which  V  is 
decomposed.  The  expanded  obstacle  A^  associated  with  A  for  a  given  orienta- 
tion 9  of  fi  is  the  pointwise  vector  difference  A  -  Bq,  where  B e  is  the  stan- 
dard placement  of  the  moving  object  B  in  which  P  lies  at  the  origin,  rotated 
by  9.  Ae  is  also  a  convex  polygonal  region  whose  sides  are  vector  (Min- 
kowski) differences  of  the  form  W  -  S  where  {W,S)  is  an  contact  pair  of 
type  I  or  n,  and  whose  vertices  have  a  similar  representation  for  contact  pairs 
(W,5)  of  type  HI. 

It  follows  from  the  results  of  [KS]  that  the  restricted  free  configuration 
space  FPq,  which  is  the  space  of  all  free  placements  of  B  having  orientation 
9 ,  can  be  represented  as  the  complement  of  the  union 

K,  =  ^U^  (AOe*  =    U^  (Ai  -  B,) 

where  A 1,  ...  ,A^  are  the  convex  polygonal  regions  into  which  V  is  decom- 
posed. See  Fig.  2.1  for  an  illustration  of  K^  and  of  FP^.  The  boundary  of 
ATe  (and  also  of  FP^)  thus  consists  of  a  collection  of  polygonal  curves  having 
finitely  many  comers.  An  edge  on  the  boundary  of  Kq  is  a  connected  portion 
of  an  edge  of  an  expanded  obstacle  induced  by  a  type  I  or  a  type  II  contact 
pair,  and  each  vertex  of  K^  is  either  a  (convex)  comer  of  an  expanded 
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obstade  induced  by  a  type  III  (X)ntact  pair,  or  a  (non-convex)  intersection 
point  of  two  sides  of  different  expanded  obstacles,  each  induced  by  a  contact 
pair  of  type  I  or  of  type  11. 

The  non-convex  corners  of  the  boundary  of  ATe  are  the  locations  of  the 
reference  point  P  at  semi  free  placc;nents  in  which  B  (at  orientation  6)  makes 
two  distinct  obstacle  contacts  simultaneously  (in  each  of  which  either  a  comer 
of  B  touches  a  wall  edge  or  a  side  of  B  touches  a  wall  comer). 

It  is  shown  in  [KS]  (cf.  also  [KLPS],  [LS2])  that,  for  each  fixed  orienta- 
tion 9  of  5,  the  number  of  non-convex  comers  oi  K^  is  only  0{n)  (actually  it 
is  only  0{m)).  An  algorithm  for  the  calculation  of  FP^,  which  makes  use  of 
this  property,  is  presented  in  [KS];  its  time  complexity  is  0(kn  log^kn), 
which  has  later  been  improved  to  0{kn  log  kn)  in  [LS2],  using  a  different 
approach  involving  generalized  Voronoi  diagrams  (and  also  in  a  recent  paper 
[BZ]  that  uses  a  line-triangle  representation  of  the  moving  object). 

The  first  step  of  our  algorithm  is  to  use  the  results  of  [KS]  or  of  [LS2]  to 
obtain  a  decomposition  of  FP^  into  connected  components,  and  to  represent 
each  such  component  Q  by  a.  discrete  connected  graph  whose  nodes  are  the 
comers  of  Q,  each  of  which  is  given  a  discrete  combinatorial  labeling  that 
does  not  depend  continuously  on  9.  (This  differs  from  the  approach  used  in 
[LSI],  namely  to  partition  FP^  into  trapezoidal  cells,  and  to  establish  adja- 
cency of  these  cells  in  FPq,  obtaining  this  way  a  connectivity  graph  CG^ 
whose  connected  components  correspond  in  a  1-1  manner  to  the  connected 
components  of  FP^.  Although  this  approach  can  be  used  in  our  case  as  well, 
it  creates  technical  difficulties  when  we  attempt  to  extend  these  cells  into  3-D 
cells  as  we  add  the  degree  of  freedom  of  rotation.  We  therefore  prefer  to  use 
the  approach  mentioned  above,  which  is  m_ore  similar  to  that  of  [SiS].) 

More  precisely,  let  Cq  =  Co(9)  be  the  collection  of  all  convex  and  non- 
convex  comers  of  K^  (and  also  of  FPe).  Define  a  vertically  maximal  corner 
of  FPq  to  be  a  convex  comer  m  of  A^e  which  has  the  largest  >'-coordinate 
among  all  points  of  AT  e  in  a  sufficiently  small  neighborhood  of  u.  For  each 
such  vertically  maximal  comer  u,  we  introduce  an  auxiliary  corner  u  * ,  which 
is  the  unique  point  on  the  boundary  of  FPo  that  lies  directly  above  u  and  is 
such  that  the  open  segment  uu*  is  wholly  contained  in  FPq.  (Note  that  since 
V  is  assumed  to  be  bounded,  FP^  is  also  bounded,  and  u  *  thus  always  exists 
and  is  well-defined.)  Let  C*  =  C*(9)  denote  the  set  of  all  such  auxiliary 
comers. 

Next  define  a  vertex  graph  VG^,  whose  set  of  nodes  is 
C  =  C(9)  =  Co  U  C*,  and  whose  edges  either  connect  pairs  of  adjacent 
comers  (including  auxihary  ones)  along  edges  of  FPq,  or  connect  vertically 
maximal  comers  u  to  their  associated  auxiliary  comers  u*  €  C*. 
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For  certain  critical  values  of  0,  the  boundary  of  FP^  may  contain  a 
comer  u  for  which  the  intersection  of  FP^  with  any  sufficiently  small  neigh- 
borhood of  u  is  disconnected;  this  would  be  the  case  if  at  orientation  6  a 
comer  of  B  touches  a  wall  comer,  while  B  also  makes  another  contact  with 
the  walls,  which  corresponds  to  a  convex  comer  of  K^  also  lying  on  another 
edge  of  Kq.  (Note  that  more  degenerate  doubk-  contacts,  in  which,  say,  two 
convex  corners  of  K^  touch  one  another,  are  ruled  out  by  our  assumptions 
that  the  obstacles  are  in  general  position;  cf.  also  [LS3].)  In  these  cases  we 
split  u  into  two  distinct  comers,  one  for  each  connected  component  of 
FPq  D  N,  for  an  arbitrarily  small  circular  neighborhood  N  of  u,  and  then  con- 
nect each  of  these  split  comers  u'  to  its  two  adjacent  corners  along  the  two 
rays  bounding  the  component  Q'  of  FPq  associated  with  u'  (cf.  Fig.  2.2). 
The  split  comer  u'  will  be  regarded  as  lying  only  on  the  boundary  of  Q' ,  and 
not  on  the  boundary  of  the  other  components  of  FPq  which  lie  near  u.  (Note 
that,  assuming  general  position  of  the  obstacles,  such  a  comer  u  cannot  be 
vertically  maximal.) 

In  the  following  subsection  we  will  consider  general  motion  of  B  includ- 
ing rotation.  It  is  important  for  analysis  of  such  general  motion  that  VGq 
does  not  depend  continuously  on  6,  but  rather  remains  invariant  except  at 
certain  finitely  many  critical  orientations,  to  be  defined  below.  For  this  reason 
we  assign  discrete  combinatorial  labelings  to  each  node  (and  thus  also  to  each 
edge)  of  VGe  in  the  following  simple  manner:  Each  convex  comer  of  A^e  is 
labeled  by  the  type  III  contact  pair  that  induces  it;  each  non-convex  corner  is 
labeled  by  the  two  contact  pairs  (of  type  I  or  II)  that  induce  it;  and  each  auxi- 
liary comer  w  *  is  labeled  as  an  auxiliary  comer  associated  with  the  type  III 
contact  pair  that  induces  the  vertically  maximal  comer  u  to  which  u* 
corresponds.  It  is  easily  checked  that  these  labelings  uniquely  define  the 
corresponding  corners  of  FPq  (at  a  given  orientation  9).  This  labeling  scheme 
turns  VGe  into  a  discrete  structure,  clearly  not  depending  continuously  on  9. 

Theorem  2.1:  Two  corners  «,  v  €  C  belong  to  the  same  connected  com- 
ponent of  VGq  if  and  only  if  they  lie  on  the  boundary  of  the  same  connected 
component  of  FPq. 

Proof:  If  the  edge  («,v)  belongs  to  VGe  then  by  definition  u  and  v  lie  on  the 
(boundary  of  the)  same  connected  component  of  FPq.  Thus  the  "only  if"  part 
of  the  Theorem  follows  by  transitive  closure.  As  to  the  "if"  part,  it  is  clearly 
true  in  the  case  in  which  u  and  v  lie  on  the  same  connected  component  of  the 
boundary  of  a  connected  component  Q  of  FP^.  U  Q  has  more  than  one  boun- 
dary component,  then  we  let  A„  (resp.  Ay)  denote  the  set  of  all  nodes  in  VGe 
reachable  from  u  (resp.  from  v)  by  a  path  in  VGe-  We  claim  that  both  A„  and 
Ay  contain  corners  lying  on  the  (unique)  exterior  boundary  E  of  Q,  so  that 
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the  preceding  argument  implies  that  A„  =  Ay,  or,  in  other  words,  u  and  v  lie 
in  the  same  cx)nnected  comporent  of  VGq.  Indeed,  suppose  that  A„  contains 
no  comer  of  E.  Let  z  €  A„  be  the  corner  with  highest  y  coordinate.  It  is  then 
easy  to  show  that  z  is  a  ver*ic2lly  maximal  corner.  (This  is  because  z  lies  on 
an  interior  boundary  of  Q,  and  if  the  upward-directed  vertical  ray  from  z  did 
not  contain  free  positions  sufficiently  near  z,  then  it  would  have  to  intersect 
the  boundary  component  containing  z  at  a  point  having  larger  y-coordinate 
than  z,  contradicting  our  choice  of  z.  This  argument  is  easy  in  2-space;  a  more 
general  argument  for  higher-dimensional  space  is  given  in  [SiS,  Section  3].) 
But  then  VGq  contains  the  edge  {z,z*),  and  z*  has  larger  y-coordinate  than  z, 
a  contradiction  which  completes  the  proof  of  the  Theorem.  □ 

2.2.  The  case  of  a  general  moilow  of  B. 

We  now  turn  to  the  general  case  in  which  B  can  both  translate  and 
rotate.  To  do  so  we  first  consider  how  the  combinatorial  characterization  of 
FPq,  provided  by  the  graph  VGq,  changes  as  9  varies. 

Definition  2.3:  An  orientation  6  of  fi  is  called  a  critical  orientation,  if  one  of 
the  following  conditions  occurs: 

(i)  There  exists  a  semi-free  placement  of  B  at  orientation  6  at  which  either  it 
makes  simultaneously  three  distinct  obstacle  contacts  involving  contact 
pairs  of  types  I  or  11,  or  it  makes  simultaneously  two  obstacle  contacts, 
one  involving  a  contact  pair  of  type  III  and  another  involving  a  contact 
pair  of  type  I  or  II.  In  other  words,  either  three  edges  of  expanded  obs- 
tacles meet  at  the  same  non-convex  corner  of  K^,  or  a.  convex  comer  of 
ATe  meets  another  edge  of  Kq  (see  Fig.  2.3(a,b)). 

(ii)  There  exists  a  vertically  maximal  convex  comer  w  of  ^e  whose  associated 
auxiliary  corner  u*  coincides  with  a  (convex  or  non-convex)  comer  of 
K,  (see  Fig.  2.3(c)). 

(iii)  (a)  Two  adjacent  edges  of  K^  become  collinear  (Fig.  2.3(d)),  or 

(b)  An  edge  of  Kq  becomes  horizontal  (see  Fig.  2.3(e)). 

Lemma  2,5:  The  vertex  graph  VGq  does  not  change  as  8  varies  in  a  suffi- 
ciently small  neighborhood  of  any  non-critical  orientation.  Furthermore,  for 
each  such  sufficiently  small  neighborhood,  the  connected  component  Q{A,Q) 
of  FPq,  corresponding  to  a  fixed  connected  component  A  of  VG^,  varies  con- 
tinuously (in  the  Hausdorff  topology  of  sets)  with  9 . 

Proof:  Filial  observe  that  as  long  as  condition  (iii)  does  not  arise,  each  edge 
or  a  convex  corner  of  an  expanded  obstacle,  corresponding  to  some  contact 
pair  O,  continues  to  appear  on  the  boundary  of  that  obstacle,  and  varies  con- 
tinuously with  9.  Moreover,  the  intersection  between  any  two  such  edges 
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(belonging  to  different  expanded  obstacles)  is  transversal,  unless  either  these 
edges  become  collinear,  or  they  meet  at  an  endpoint  of  one  of  them.  Hence 
as  long  as  conditions  (i)  and  (iii)  do  not  arise,  each  such  intersection  is 
transversal,  and  thus  also  varies  continuously  with  G.  Also,  as  long  as  no 
three  edges  of  expanded  obstacles  meet  at  the  same  non-convex  comer  of  ATe , 
and  no  convex  comer  of  Kq  meets  an  edge  of  that  set,  it  follows  that  as  we 
vary  9,  each  non-convex  comer  of  K^,  formed  by  intersection  of  edges 
corresponding  to  two  contact  pairs  Oi,  O2,  will  continue  to  appear  on  the 
boundary  of  K^,  continue  to  be  induced  by  the  same  two  contact  pairs,  and 
vary  continuously  with  0.  It  follows  that  as  6  varies  slightly,  no  connected 
subcomponent  of  FPq  (or  even  no  component  of  the  boundary  of  FPq) 
shrinks  to  a  point  and  disappears,  no  such  component  newly  appears,  no  two 
such  components  merge  into  a  single  component,  nor  does  such  a  component 
split  into  two  subcomponents.  Furthermore,  each  connected  component  of 
the  boundary  of  FPq  retains  the  same  combinatorial  representation  as  a  circu- 
lar sequence  of  convex  and  non-convex  comers  induced  by  the  same  combina- 
tions of  obstacle  contacts,  and  each  of  these  comers  varies  continuously  with 
9.  This  also  implies  that  each  connected  component  of  FPq  varies  continu- 
ously with  0. 

To  complete  the  proof  of  the  Lemma,  let  u  be  a  vertically  maximal 
comer  of  Kq.  First  note  that  as  long  as  we  do  not  cross  any  critical  orienta- 
tion, u  remains  vertically  maximal  as  9  varies.  Also,  as  long  as  condition  (ii) 
does  not  arise,  the  auxiliary  comer  u  *  associated  with  u  remains  on  the  same 
edge  of  ATe  between  the  same  two  vertices  of  VG^,  and  varies  there  continu- 
ously with  e.  All  these  observations  clearly  imply  that  VGq  remains  constant 
as  9  varies  through  non-critical  orientations.   □ 

Critical  orientations  of  type  (i)  were  analyzed  in  [LS3],  where  it  was 
shown  that  there  are  at  most  O (knX(,{kf^))  such  orientations.  Section  3 
presents  an  algorithm  for  the  calculation  of  these  orientations,  which  runs  in 
time  0{kn\^(kn)  log  kn). 

Critical  orientations  of  type  (iii)  are  trivial  to  analyze.  In  fact,  there  are 
only  0{kn)  such  orientations,  whose  calculation  is  straightforward. 

Finally,  critical  orientations  of  type  (ii)  require  involved  analysis,  similar 
to,  but  somewhat  simpler  than  that  in  [LS3].  This  analysis  will  be  presented 
in  Section  3  below,  where  we  show  that  the  number  of  such  critical  orienta- 
tions is  also  at  most  OikriX^ikn)),  and  that  they  can  also  be  calculated  in  time 
0(knX(,(kn)  log  kn). 

Suppose  then  that  we  have  calculated  all  these  critical  orientations,  and 
sorted  them  in  circular  order.  Let  T  denote  the  sorted  collection  of  all  such 
orientations. 
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Let  00  be  a  noncritical  orientation,  and  let  u  be  (the  discrete  labeling  of) 
a  comer  in  VGqq.  We  will  associate  with  u  a  life-span  L„  =  (61,62),  which  is 
the  maximal  (open)  interval  containing  60  such  that  u  ^  VGo  for  all 
6  ^  (61,62),  and  such  that  u  does  not  become  coincident  with  any  other 
comer  of  VGe  for  any  such  6.  We  also  associate  a  similar  life-span  interval 
L^  CLuHLy  with  each  edge  e  =  («,v)  of  VGqq,  defined  in  a  similar 
manner.  We  call  the  pair  (m,L„)  an  extended  comer,  and  the  pair  (e,Lj  an 
extended  edge. 

We  next  construct  an  (extended)  edge  graph  EG,  whose  nodes  are  all 
these  extended  comers,  and  each  of  whose  edges  connects  an  extended  corner 
(m,L„)  to  another  extended  corner  (v,Lv),  if  (m,v)  is  an  edge  of  VGe  for  some 
e€L„nLv  (i.e.  if  the  life  span  of  the  edge  (u,v)  is  non-empty). 

The  graph  EG  is  used  to  represent  FP  as  follows.  For  each  (label  of  a) 
comer  of  some  Kq,  the  node  (u,L„)  of  EG  represents  an  edge  of  the  boun- 
dary of  FP  which  consists  of  all  placements  Z  =  (X,  6)  of  B  such  that  6  €  L„ 
and  such  that  at  placement  Z  either  B  makes  a  double  contact  involving  the 
two  contact  pairs  inducing  u,  if  u  is  a  non-convex  corner  of  K^,  or  B  makes 
an  obstacle  contact  a  contact  involving  the  type  III  contact  pair  inducing  u,  if 
w  is  a  convex  comer  of  ATe;  if  «  is  an  auxiliary  comer,  then  {u,Lu)  represents 
an  arc  along  a  face  of  FP,  consisting  of  placements  {X,  6)  at  which  B  makes  a 
single  obstacle  contact  induced  by  some  contact  pair  of  type  I  or  n,  such  that 
X  lies  directly  above  the  convex  comer  of  ATe  inducing  u. 

It  is  also  worth  noting  that  by  definition  each  edge  e  of  FP  is  the  inter- 
section of  two  surfaces  bounding  FP,  and  is  thus  the  locus  of  placements  Z  of 
B  at  which  it  satisfies  two  independent  constraints  involving  contacts  with 
obstacles.  Hence  at  these  placements  either  B  makes  two  distinct  contacts 
with  the  walls,  involving  two  distinct  contact  pairs  Oi  =  {Wi,Si), 
02  =  iW2,S2)  of  type  I  or  II,  or  it  makes  a  single  contact  of  one  of  its  comers 
S  against  a  wall  comer  W,  i.e.  involving  the  type  III  pair  O  =  {W,S).  It  fol- 
lows that  e  must  correspond  to  some  edge  of  EG  defined  either  in  terms  of 
the  non-convex  comer  labeled  by  Oi,  O2,  or  by  the  convex  corner  labeled  by 
O,  except  in  degenerate  situations  in  which  all  placements  along  e  have  the 
same  orientation  6;  these  can  arise  in  certain  degenerate  situations,  including 
e.g.  the  following  cases  (cf.  Proposition  2.1  of  [LS3]): 

(a)  Wi,  W2  are  both  wall  comers  and  ^i  =52;  (see  Fig.  2.4  (a)). 

(b)  Wi,  W2  are  both  wall  comers  and  Si,  S2  are  parallel  sides  of  B; 

(c)  Wi,  W2  are  parallel  wall  edges;  (see  Fig.  2.4  (b)). 

(d)  Wi  is  a  wall  edge  and  52  is  a  side  of  B  parallel  to  Wi . 
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Nevertheless,  we  will  show  that  the  edges  of  FP  that  are  represented  as  nodes 
of  EG  are  sufficient  for  adequate  represantation  of  FP. 

The  edges  of  EG  also  have  natural  interpretation  in  terms  of  FP.  Let 
i  =  [(m,L„),  iy,Ly)]  be  an  edge  of  EG,  and  let  e.^,  e^  denote  the  two  edges  of 
FP  corresponding  to  (w,L„)  and  (v.L^)  respex:tivdy.  Then  ^  represents  adja- 
cency of  e„  and  e^  along  a  face /of  FP;  more  precisely,  for  each  6  C  ^(^.v). 
the  cross-section  of  /  at  9  contains  a  straight  segment  whose  endpoints  are 
(labeled  by)  u  and  v.  An  exception  is  when  v  =  u*;  then  ^  represents  in  a 
similar  manner  adjacency  of  e^  and  e^*  along  an  auxiliary  "vertical"  surface 
inside  FP. 

We  claim  that  EG  captures  the  connectivity  of  FP  in  the  following  sense: 

Proposition  2.6:  (A)  To  each  connected  component  A  of  EG  there 
corresponds  a  unique  connected  component  Q  of  FP,  such  that  the  nodes  in  A 
represent  precisely  those  edges  of  FP  which  bound  Q  (and  which  do  not  lie  in 
a  single  6 -cross  section  of  FP). 

(B)  Each  connected  component  of  FP  has  at  least  one  bounding  edge  that  is 
represented  by  a  node  of  EG. 

Proof:  The  proof  of  (B)  is  straightforward.  Indeed,  let  (2  be  a  connected 
component  of  FP,  and  let  Z  =  (AT, 9)  ^  Q.  The  cross  section  (2e  =  G  H  FP^ 
of  Q  contains  X  and  has  thus  nonempty  interior.  Since  V  is  assumed  to  be 
bounded,  so  must  be  j2e-  Thus  j2e  has  a  bounded  polygonal  boundary,  which 
therefore  must  contain  at  least  one  non-convex  comer  u  of  K^  (e.g.  choose  u 
to  be  the  corner  of  Ce  with  the  largest  y  coordinate).  But  then  it  is  clear  that 
(m,L„)  is  a  node  of  EG  which  represents  an  edge  of  Q. 

The  proof  of  (A)  consists  of  two  parts.  In  the  first  part  we  argue  that  if 
(u,L„),  {v,Ly)  are  two  adjacent  nodes  of  EG,  then  they  represent  edges  of  FP 
which  bound  the  same  connected  component.  This  claim  follows  immediately 
from  definition  of  EG.  Thus  it  follows  by  transitive  closure  that  all  the  nodes 
in  a  connected  component  of  EG  represent  edges  of  FP  bounding  the  same 
connected  component  of  FP. 

The  second  part  proves  the  converse  statement,  namely  that  any  pair  of 
edges  of  FP  which  bound  the  same  connected  component  of  FP,  and  neither 
of  which  is  contained  in  a  single  9-cross  section  of  FP,  are  represented  by 
nodes  of  EG  which  belong  to  the  same  connected  component  of  this  graph. 
The  proof  of  this  statement  consists  of  the  following  steps. 

The  set  of  all  critical  orientations  partitions  the  angular  space  into  dis- 
joint open  non-critical  intervals.  By  our  assumptions  that  B  and  the  obstacles 
are  in  general  position,  it  follows  that  at  each  critical  orientation  only  a  single 
criticality  occurs.  Define  a  decomposition  of  FP  into  disjoint  connected  cells 
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as  follows.  Let  /  be  a  non-critical  interval,  and  let  A  be  a  connected  com- 
ponent of  VGe  for  any  (hence  all)  9  €  /.  The  cell  c  =  c{I,A)  of  FP  associ- 
ated with  /  and  A  is  defined  as 

c(/,A)  =  {(X,e):G  €/,     XC(2(A,e)}, 

where,  as  above,  QiA,Q)  is  the  (unique)  connected  comjx)nent  of  FP©  whose 
boundary  comers  appear  in  A. 

It  is  clear  that  the  cells  defined  in  this  manner  are  open,  connected,  and 
pairwise  disjoint,  and  that  the  union  of  their  closures  covers  the  entire  space 
FP. 

Two  cells  c  =  c{I,A),  c'  =  c(/',A')  are  called  adjacent  if  the  intervals  /, 
/'  have  a  common  endpoint  6*  and  there  exists  a  free  placement 
{X*  ,Q*)  ^  FP  lying  in  the  closures  of  both  cells  c,c' . 

Now  let  Z  ={X,  9),  Z  =  (A", 9')  be  two  free  placements  lying  in  the  same 
connected  component  of  FP,  and  let  p  {t)  be  a  continuous  path  in  FP  connect- 
ing these  two  placements.  Without  loss  of  generality  we  may  assume  that 
p{t)  crosses  between  the  cells  of  FP  as  defined  above  only  finitely  often, 
intersecting  each  cell  boundary  transversally,  and  in  fact  enters  each  such  cell 
at  most  once.  Lei  c  i ,  .  .  .  ,  c,  be  the  sequence  of  cells  traversed  by  p  {t)  in 
this  order.  It  is  then  clear  that  each  pair  of  cells  c/,  c/  +  i  in  this  sequence  are 
adjacent  in  the  sense  just  defined.  Conversely,  if  c,  c'  are  adjacent  cells  then 
one  can  show  that  there  exists  a  free  motion  of  B  from  any  placement  in  c  to 
any  placement  in  c',  which  is  contained  in  the  union  of  the  closures  of  these 
two  cells  and  which  intersects  the  boundary  between  c  and  c'  transversally  at 
a  single  point.  (To  show  this,  one  has  to  move  B  from  a  placement  in  c  (or 
in  c')  to  a  free  placement  (X*,9*)  lying  in  both  closures  of  c  and  of  c';  to 
argue  that  such  a  motion  is  always  possible  without  having  to  enter  other 
cells,  we  use  a  straightforward  generalization  of  the  argument  given  in  [LSI, 
Section  2.2];  in  fact,  as  argued  there,  the  required  motion  can  be  achieved  by 
a  pure  rotation  around  (X*,9*).) 

We  can  therefore  define  a  connectivity  graph  CG,  whose  nodes  are  the 
cells  in  the  above  decomposition  of  FP,  and  whose  edges  connect  pairs  of 
adjacent  cells.  The  above  argument  then  implies  that  two  cells  c,  c'  lie  in  the 
same  connected  component  of  CG  if  and  only  if  they  are  contained  in  the 
same  connected  component  of  FP. 

To  establish  our  claim  it  is  thus  sufficient  to  prove  that 

(a)  The  boundary  of  each  cell  c  of  FP  contains  at  least  one  edge  which  is 
contained  in  an  edge  e^  represented  by  some  node  (w,L„)  of  EG. 

(b)  If  two  cells  c  =  c{I,A),  c'  =  c(/',A')  of  FP  are  adjacent,  and  if  (u,L„), 
(v,L»,)  are  two  nodes  of  EG  whose  corresponding  edges  e^,  e^  intersect 


-13- 

the  boundaries  of  c,  c'  respectively,  then  they  lie  in  the  same  connected 
component  of  EG.  Moreover,  the  same  propert}^  also  holds  if  c  =  c'. 

To  prove  property  (a),  let  c  =  c  {I, A)  be  such  a  cell  of  FP,  and  let  m  be  a 
node  in  A.  Then,  for  any  Q  i  I,  u  represents  a  comer  of  the  connected  com- 
ponent G  =  (2  (A,  e)  of  F/*e .  It  thus  follows  by  definition  that  /  is  a  subset  of 
L„  and  that  e„  has  the  desired  property. 

As  to  property  (b),  we  first  claim  that,  given  two  adjacent  cells 
c=c(/,A),  c'  =  c(/',A'),  there  exists  a  node  {w,L^)  of  EG  representing  an 
edge  of  FP  which  meets  the  boundaries  of  both  c  and  c' .  Indeed,  let  0*  be 
the  common  endpoint  of  /  and  /'.  By  assumption,  6*  involves  only  one  criti- 
cality.  If  this  criticality  does  not  affect  the  components  A,  A'  of  VG^,  then 
clearly  they  must  be  equal,  and  any  comer  w  in  A=A'  yields  a  node  {w,L^) 
of  EG  with  the  desired  property. 

Suppose  next  that  at  9*  a  criticality  which  affects  A  or  A'  occurs.  Let 
G  =  (2(A,e)  fore  €  /(resp.  Q'  =  Q{A' ,B')  for  Q'  ^  /')  be  the  connected  com- 
ponent of  FPe  (resp.  of  FP^')  represented  by  A  (resp.  by  A').  Note  that  both 
Q  and  Q'  are  bounded  polygonal  regions,  so  that  the  exterior  boundary  of 
each  must  have  at  least  three  corners.  Furthermore,  each  criticality  can  affect 
at  most  two  corners  of  Q  (resp.  of  Q').  Indeed,  the  only  type  of  criticality 
that  might  affect  three  comers  of,  say,  Q  is  that  in  which  three  edges  of  Kq 
meet  at  a  common  point  at  9  =  6*.  But  this  can  affect  three  comers  of  Q  only 
if  2  is  bounded  by  just  this  three  edges,  and  thus  shrinks  to  a  point  as  9  -  0* ; 
this  however  is  impossible  because  then  c  and  c'  could  not  have  been  adjacent 
in  FP. 

These  arguments  are  easily  seen  to  imply  that  A  and  A'  must  share  a 
common  node  w  which  is  not  affected  by  the  critical  change  occurring  at  0  * , 
and  it  then  follows  by  definition  that  (w,L„)  meets  the  boundaries  of  both  c 
and  c' ,  as  asserted.  It  is  thus  sufficient  to  prove  the  second  claim  in  (b), 
because  the  first  claim  will  then  follow  by  combining  the  facts  that  (u,L„)  and 
{w,L„)  (as  edges  bounding  c)  lie  in  the  same  connected  component  of  EG, 
and  that  (v,Lv,)  and  {w,L^)  (as  edges  bounding  c')  also  lie  in  the  same  con- 
nected component  of  EG. 

Thus  let  (u,L„),  {w,L^)  be  two  nodes  in  EG  representing  edges  of  FP 
both  meeting  the  boundary  of  the  same  cell  c  =  c{I,A).  Then  it  follows 
from  Theorem  2.1  that  u  and  w  belong  to  the  same  connected  component  of 
VGe,  for  each  9  6  /,  so  that,  by  definition,  (w,L„)  and  {w,L^)  belong  to  the 
same  connected  component  of  EG.  The  proof  of  (b),  and  hence  also  of  part 
(A)  of  our  proposition,  is  thus  completed.  □ 

Defmition  2.4:  Using  the  observations  made  in  the  proof  of  (B)  above, 
we  can  also  define  a  map  <t)(Z)  which  maps  each  Z  ^  FP  to  a  node  of  EG 
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representing  an  edge  of  FP  which  bounds  the  connected  component  of  FP 
containing  Z.  specifically,  given  any  free  placement  Z  =  (X,6)  of  B,  we  can 
translate  B  upwards  from  Z  {if:,  in  the  positive  y  direction)  imtil  it  reaches  a 
new  placement  Z*  =  (X*,9>  at  which  it  makes  contact  with  an  obstacle.  If 
this  contact  involves  a  contact  pair  of  type  HI,  or  consists  of  two  simultane- 
ous contacts  involving  pairs  of  types  I  or  H,  then  AT*  is  a  comer  of  FP^,  and 
we  let  <t)(Z)  be  the  node  of  EG  induced  by  that  comer.  Otherwise,  we  con- 
tinue to  translate  B  from  Z*  leftwards,  maintaining  the  obstacle  contact  that 
involves  the  same  contact  pair,  until  we  reach  a  placement  Z**  =  (X**,0) 
for  which  X**  is  a  corner  of  FPq,  and  then  continue  as  above.  Again,  the 
fact  that  FPq  is  bounded  implies  that  Z*,  Z** ,  and  thus  also  ct)(Z),  are 
always  well  defined. 

We  thus  obtain  the  following  main  theorem  of  this  section: 

Theorem  2.2:  Let  Z  ,  Z'  €FP  be  two  free  placements  of  B.  Then  there  exists 
a  collision-free  motion  of  B  from  Z  to  Z'  if  and  only  if  the  nodes  <^{Z),  ^{Z') 
of  EG  as  defined  above  belong  to  the  same  connected  component  of  EG. 

3.  Algorithmic  Details 

We  now  turn  to  describe  an  efficient  implementation  of  our  motion  plan- 
ning solution.  We  first  present  an  algorithm  that  computes  VG^  for  any  fixed 
6;  this  is  an  expansion  of  the  algorithm  described  in  [KS]  for  the  calculation 
of  FP  in  the  case  of  purely  translational  motion  of  a  convex  polygonal  body 
amidst  convex  polygonal  obstacles.  In  the  next  Subsection  we  show,  in  a 
manner  similar  to  [LS3],  that  the  number  of  critical  orientations  of  type  (ii)  is 
0(kn  \^{kn)).  We  follow  with  a  detailed  description  of  the  algorithms  that 
find  all  critical  orientations,  and  then  construct  the  extended  edge  graph  EG. 
Finally,  using  these  structures,  we  show  how  to  plan  continuous  obstacle- 
avoiding  motion  of  B  between  any  specific  initial  and  final  placements. 

3.1.   Constructing  VG^ 

Let  e  be  a  non-critical  orientation.  The  following  algorithm  constructs 
VGq  in  four  stages.  The  first  three  stages  produce  FP^  in  essentially  the  same 
manner  as  in  [KS],  and  the  additional  fourth  stage  computes  VGq  from  FPo. 

(a)  For  each  of  the  convex  polygonal  obstacles  A  into  which  V*^  is  decom- 
posed (see  Subsection  2.1),  calculate  the  expanded  obstacle  Aq  =  A  -  Bq  in 
time  linear  in  the  sum  of  the  number  of  corners  of  A  and  of  B  (using  e.g  the 
technique  of  [GRS];  note  that  the  number  of  comers  of  Ae  is  also  bounded 
by  that  sum  [GRS]).  Thus  the  total  number  of  convex  comers  of  all  the 
expanded  obstacles  is  0{km+n)  and  they  can  all  be  calculated  in  time 
0{km  +n),  where  m  is  the  number  of  the  convex  subparts   of  V,  and  k,n  art 
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as  above.  Since  m<n  we  can  denote  this  complexity  as  0{kn).  However,  if 
m«n  then  0{kn)  may  be  a  gross  overestimation  of  the  actual  complexity  of 
the  expanded  obstacles  at  any  given  8.  See  below  for  a  discussion  on  how  this 
observation  affects  the  overall  complexity  of  our  algorithm. 

(b)  Calculate  the  boundary  of  ATe,  which  is  the  union  of  all  these 
expanded  obstacles.  To  do  this  we  use  a  technique  due  to  [OWW]  which  cal- 
culates the  boundary  of  the  union  of  several  superimposed  polygonal  planar 
regions  in  time  0((m  +r)  log  n),  where  n  is  the  number  of  line  segments  con- 
stituting the  boundaries  of  these  regions,  and  /  is  the  number  of  intersections 
between  these  segments.  We  combine  this  technique  with  the  following 
divide  and  conquer  approach: 

1.  Recursively  find 

lm/2] 

G  =    ^U    (AOe* 

m 

//  =         U        (Ai)t 

1=  [mil  +ll 

2.  Find  the  contour  of  Kq  =  GUH  using  the  technique  of  [OWW]. 

Since  each  point  of  intersection  of  the  boundaries  of  G  and  H  must  lie  on  the 
boundary  of  their  union,  it  follows  by  a  theorem  proved  in  [KS],  that  the 
number  of  these  intersection  points  is  only  0{m),  so  that  step  2  of  the  above 
algorithm  runs  in  time  0{kn  log  kn),  and  thus  the  entire  algorithm  runs  in 
time  0{kn  log  kn  log  m). 

(c)  Find  the  connected  components  of  FP^  and  represent  them  as  an 
"inclusion  tree",  each  of  whose  nodes  is  a  connected  component  of  the  boun- 
dary of  FPq,  as  follows  [KS].  Note  that  FP^  is  a  general,  possibly  discon- 
nected polygonal  planar  region,  whose  components  need  not  be  simply  con- 
nected. However,  since  0  is  assumed  to  be  non-critical,  each  connected  com- 
ponent of  the  boundary  of  FP^  must  be  a  closed  simple  polygonal  path.  In 
the  inclusion  tree  that  we  construct  each  node  represents  a  simple  connected 
component  of  this  boundary,  and  such  a  component  C2  is  a  child  of  another 
Ci  if  C2  is  contained  in  Ci,  and  the  region  between  them  is  either  contained 
entirely  within  FP^  or  is  disjoint  from  FPq.  The  root  of  the  inclusion  tree  is 
a  nominal  "infinite  boundary"  consisting  of  all  points  at  infinity.  Each  con- 
nected component  C  of  FP^  is  represented  in  the  inclusion  tree  by  the  node 
which  represents  the  exterior  polygonal  boundary  of  C,  and  by  all  children  of 
that  node  (if  any),  which  represent  the  connected  components  of  the  interior 
boundary  of  C.  The  entire  tree  is  easily  calculated  in  time  0(kn  \og  kn), 
using  a  straightforward  sweeping  technique  (cf.  [KS]  for  more  details). 
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(d)  Finally  we  construct  VGe  from  FP©  as  follows.  Initially,  we  take  all 
comers  of  FPe  to  be  vertices  of  VGe,  and  connect  them  in  VG^  along  edges 
of  FP^.  Then,  sweeping  a  vertical  scan  line  across  the  plane,  we  can  easily 
calculate  the  auxiliary  comer  associatexi  with  each  vertically  maximal  corner 
of  ATe,  in  overall  time  0{kn  log  kr).  We  then  add  each  of  these  auxiliary 
comers  m*  as  a  vertex  of  VGe,  split  the  edge  e  of  FP©  containing  u*  into 
two  subsegments  at  w  * ,  remove  e  from  VG  e ,  and  finally  connect  w  *  to  the 
vertically  maximal  comer  u  that  induced  it,  and  to  the  two  endpoints  of  e. 

Clearly  the  time  complexity  of  the  entire  algorithm  is 
0{kn  log  kn  log  m). 

Lemma  3.1:  The  vertex  graph  VGq  is  planar,  has  0{kn)  edges  and  comers 
and  each  vertex  in  VGe  is  incident  to  at  most  three  edges. 

Proof:  The  planarity  follows  immediately  from  the  fact  that  FPq  is  planar 
and  the  added  vertical  segments  which  emanate  upwards  from  vertically  max- 
imal comers  do  not  intersect  any  other  segment  except  at  a  vertex  of  VG  e . 

The  number  of  corners  and  edges  of  FPq  is  0{kn)  and  the  number  of  auxih- 
ary  comers  does  not  exceed  the  number  0{kn)  of  convex  corners  of  FFe,  and 
each  such  auxihary  corner  adds  two  edges  to  VGe-  Hence  VGe  has  0{kn) 
edges  and  comers. 

Finally,  it  is  plain  that  each  vertex  in  VGe  is  incident  to  two  edges,  except  for 
the  vertically  maximal  vertices  and  their  corresponding  auxihary  vertices, 
each  being  incident  to  three  edges.  □ 

3.2.  The  number  of  critical  orientations  of  type  (ii). 

We  next  turn  to  the  problem  of  estimating  the  number  of  critical  orienta- 
tions of  type  (ii).  Recall  that  these  orientations  arise  when  the  auxiliary 
comer  u  *  associated  with  a  vertically  maximal  comer  u  of  Kq  coincides  with 
a  convex  or  a  non-convex  comer  oi  K^.  Also  recall  that  a  vertically  maximal 
comer  u  is  induced  by  a  contact  pair  O  =  (W,^),  where  W  is  a  wall  corner 
and  S  a  comer  of  B  such  that  u  has  the  largest  ^-coordinate  among  all  points 
of  ATe  in  a  sufficiently  small  neighbourhood  of  u. 

The  number  of  critical  orientations  8  at  which  such  u  lies  immediately 
below  a  convex  comer  v  of  K^  is  0{k^n^)  and  a  superset  of  these  orienta- 
tions, also  of  size  0{k^n^),  can  be  calculated  in  a  straightforward  manner  by 
iterating  over  all  0{k^n^)  possible  two  contact  pairs  of  type  HI.  We  thus 
proceed  to  estimate  the  number  of  critical  orientations  at  which  u  lies  directly 
below  a  non-convex  comer. 
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Let  u  be  such  a  fixed  vertically  maximal  comer  induced  by  a  type  m  con- 
tact pair  O  =  (W,5).  For  each  contact  pair  Ci  =-  {Wi,Si)  of  type  I  or  n, 
and  for  each  6  define  the  function  Gqoi  (6)  '^  be  tlie  dii^tance  along  the  verti- 
cal ray  from  the  comer  u  upwards  to  its  intersection  with  the  edge  induced  by 
the  pair  Oi  if  such  an  intersection  exists,  and  +^  olhe^rwise. 

Let  n  =  Ylooi  denote  the  domain  of  definition  of  Gooii^)-   Note  that  11 
need  not  in  general  be  connected,  but,  as  in  [LS3],  we  have: 
Lemma  3.2:  ll^oj  consists  of  at  most  4  subintervals. 

Proof:  Fix  the  comer  5  of  fi  at  the  wall  comer  W  and  let  B  rotate  about  that 
common  comer.  Each  point  X  in  B  then  traces  a  circular  arc,  denoted  by  Cx , 
about  W.  Suppose  first  that  5i  is  a  corner  of  B  (snd  that  Wi  is  a  wall  edge). 
It  is  clear  that  Gqo li^)  is  defined  in  this  case  if  and  only  if  the  vertical  ray 
emanating  upwards  from  5i  intersects  Wi,  and  that  Gooi(9)  is  equal  to  the 
length  of  the  vertical  segment  from  Si  to  the  intersection  point  of  this  ray 
with  Wi.  Similar  considerations  apply  when  Si  is  a  side  of  B  (and  Wi  is  a 
wall  comer).  In  this  case  Gqo li^)  is  defined  if  and  only  if  the  vertical  ray 
emanating  downwards  from  Wi  intersects  Si,  and  as  before  Gooj(9)  is  then 
equal  to  the  distance  between  Wi  and  the  intersection  point  of  that  ray  with 
Si. 

Thus  Goo li^)  is  defined  for  orientations  0  of  5  that  fulfill  the  following 
three  conditions: 

(1)  The  corner  u  induced  by  O  exists  and  is  vertically  maximal. 

(2)  B  can  make  at  orientation  6  an  obstacle  contact  involving  the  contact  pair 
Oi  =  iWuSi). 

(3)  The  vertical  ray  emanating  upwards  from  u  intersects  the  expanded  obs- 
tacle edge  induced  by  C>i. 

Thus  the  domain  of  Goo li^)  is  the  intersection  of  the  three  domains 
1 1,  1 2,  1 3  satisfying  respectively  these  three  conditions. 

A  necessary  condition  for  Ii,  the  domain  of  orientations  for  which  u 
exists  and  is  vertically  maximal,  to  be  nonempty  is  that  W  itself  be  vertically 
maximal,  in  the  sense  that  it  is  the  point  with  the  largest  y  coordinate  among 
all  points  of  V  in  some  small  neighborhood  of  W.  Assuming  that  this  is  the 
case,  it  is  then  clear  that  the  angular  interval  where  u  exists  and  is  vertically 
maximal  must  be  contained  in  the  interval  of  orientations  of  B  whose 
corresponding  positions  of  P  on  the  circle  Cp  lie  in  the  upper  half  of  that  cir- 
cle. 

It  is  also  clear  that  the  domain  of  orientations  at  which  a  contact  induced  by 
O  is  possible,  is  an  angular  interval  whose  length  is  equal  to  the  exterior 
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angle  at  W  minus  the  angle  of  6  at  S.  Thus  /i  is  an  angular  interval  of  length 
less  than  or  equal  to  tt.   (see  Fig.  3  1(a)). 

It  is  also  clear  that  Ii,  the  domain  of  orientations  at  which  the  contact  O^ 
is  possible,  is  an  angiilar  interval  whose  length  is  equal  either  to  the  exterior 
angle  of  B  at  5i  minus  it  if  5 1  is  a  vertex  of  B,  or  to  ir  minus  the  convex 
angle  at  W^,  if  W^  is  a  wall  comer.  Thus  the  intersection  of  /i  and  1 2  is  an 
angular  interval  of  length  strictly  less  than  it  (see  Fig.  3.1(b)). 

As  to  73 ,  suppose  first  that  Wi  is  a  wall  edge.  Let  l^  (resp.  /2)  be  a  vert- 
ical ray  emanating  from  the  rightmost  (resp.  leftmost)  endpoint  of  W-^  down- 
wards, and  consider  the  infinite  trapezoidal  strip  bounded  by  Wi,  l^,  and  I2. 
It  is  easily  checked  that  !■>,  consists  of  these  orientations  of  B  at  which  Si  Hes 
within  that  trapezoid,  or,  in  other  words,  each  connected  component  of  I-^ 
corresponds  to  a  connected  component  of  the  intersection  of  Cs^  with  this 
trapezoid.  Thus  in  this  case  I2,  consists  of  at  most  three  subintervals.  (See 
Fig.  3.1(c)).  Furthermore,  if  1 2,  does  consist  of  three  subintervals,  then  each 
of  them  has  length  at  most  tt  . 

Suppose  next  that  Wi  is  a  comer,  and  let  /^  be  a  vertical  ray  emanating 
from  Wi  downwards.  Then  73  consists  of  these  orientations  of  B  at  which  the 
side  5 1  intersects  l^.  Here  two  conditions  need  be  met  simultaneously.  First, 
Wi  must  lie  above  the  hne  containing  5i  and  W,  W^  must  lie  on  different 
sides  of  this  hne,  and  second,  the  two  endpoints  of  S\  must  lie  on  different 
sides  of  /i.  As  to  the  first  condition,  it  is  easily  checked  that  it  can  be 
expressed  as  a  hnear  inequality  in  sin  6  and  cos  9  of  the  form 
p  sin  (6  — a)  >  q,  where  p  and  q  are  both  positive,  and  can  therefore  have  a 
solution  set  consisting  of  at  most  one  angular  interval  of  length  at  most  tt 
(see  Fig.  3.1(d)).  The  second  condition  can  change  from  being  tme  to  being 
false  or  vice  versa  whenever  one  of  the  endpoints  of  5i  crosses  the  line  con- 
taining /i.  Thus  there  are  at  most  four  orientations  at  which  the  second  con- 
dition can  change,  again  yielding  at  most  two  angular  intervals  at  which  this 
second  condition  can  hold.  Thus,  intersecting  the  first  interval  with  the  latter 
two,  we  conclude  that  73  consists  of  at  most  three  angular  intervals  in  this 
case  too,  where  each  of  these  subintervals  is  of  length  at  most  tt. 

It  follows  that  in  both  cases  7i  072  073  consists  of  at  most  three  angular 
intervals.  Transforming  these  intervals  to  the  range  0  <  6  <  2-77,  we  con- 
clude that  the  domain  rio^j  of  Gooi(0)  consists  of  at  most  four  subintervals 
in  this  linearized  range.   This  completes  the  proof  of  the  lemma.  □ 

If  IIooj  is  indeed  not  connected  we  will  consider  each  connected  portion 
of  ^ooi  (6)  as  a  seperate  partially  defined  function  and  will  still  have  at  most 
0{kn)  such  functions  for  each  fixed  contact  pair  O. 
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We  next  define  the  lower  envelope  of  the  set  of  functions  Goot  ^^t  O: 

^o{Q)  =  min  Goo,(e). 

Note  that  if  /  is  a  (maximal)  interval  over  wh'ch  ^I'.-i  is  attained  by  a  single 
function  Goop  then  for  each  6  6/  the  auxiliary  ccmrj  u  *  lies  on  the  expanded 
obstacle  edge  induced  by  O/.  Furthermore,  each  breakpoint  6  on  the  envelope 
^o(9)  (i.e.  a  point  at  which  ^o(6)  is  attained  simultaneously  by  two  func- 
tions Gooi,  Gooj)  is  an  orientation  in  which  the  edge  of  K^  containing  u* 
changes  (from  the  one  labeled  by  0/  to  that  labeled  by  Oj);  conversely,  every 
such  change  in  the  edge  containing  u*  corresponds  to  a  breakpoint  of  '^q{Q). 

Next  we  show  that  each  pair  of  the  functions  Gooi  intersect  in  at  most 
four  points: 

Lemma  3.3:  Let  0  =  (W,S)  be  a  contact  pair  inducing  a  vertically  maximal 
comer  u  and  let  O^,  O2  be  two  other  distinct  contacts  pairs.  Then  there  are 
at  most  four  orientations  0  in  which  u  lies  directly  below  the  non-convex 
comer  formed  by  the  intersection  of  the  expanded  obstacle  edges  induced  by 
Oi  and  O2  respectively. 

Proof:  Follows  immediately  from  the  analysis  in  [SS]  and  Lemma  2.2  in 
[LS3].  Specifically,  consider  the  curve  701O2  traced  by  5  as  5  makes  simul- 
taneously two  obstacle  contacts  involving  the  pairs  O^  and  O2.  As  shown  in 
[SS],  this  curve  is  either  a  straight  segment  or  part  of  an  ellipse  or  part  of  a 
quartic  algebraic  curve.  The  desired  placement  at  which  u  lies  directly  below 
the  intersection  point  of  the  edges  induced  by  O^  and  02  is  the  intersection 
between  this  curve  and  the  vertical  line  emanating  upwards  from  W.  Since 
yoi02  i^  ^t  most  quartic  there  are  at  most  four  such  intersections.  □ 

It  therefore  follows  (cf .  [At])  that  the  number  of  breakpoints  along  ^  is 
at  most  X6(^")-  To  see  this,  it  suffices  to  show  that  each  pair  gi  =  Gooi, 
8j  —  Goo  of  these  partial  functions  can  alternate  along  the  lower  envelope 
^o  ^t  most  7  times.  Indeed,  the  two  extreme  portions  in  such  an  alternation 
can  occur  within  the  symmetric  difference  of  the  domains  of  gi  and  of  gj 
(observe  that  this  symmetric  difference  consists  of  at  most  two  intervals,  on 
each  of  which  just  one  of  these  functions  is  defined),  and  at  most  five  addi- 
tional altemating  portions  of  the  graphs  of  gi,  gj  can  occur  along  the  lower 
emvelope  within  the  intersection  of  the  domains  of  these  two  functions, 
because  gi  and  gj  must  intersect  between  each  pair  of  such  adjacent  altemat- 
ing occurrences,  and  by  Lemma  3.3  only  four  such  intersections  can  occur. 
Since  the  number  of  vertically  maximal  comers  is  at  most  0{kn),  we  have  in 
summary 
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Theorem  3.1:  The  total  number  of  critical  orientations  of  type  (ii)  is  at  most 

0{kn\6ikn)). 

An  efficient  procedure  for  calculating  ihese  orientations  will  be  described  in 

the  following  Subsections. 

3.3.   Calculation  of  all  critical  orientations. 

We  now  proceed  to  describe  a  collection  of  algorithms  that  actually  cal- 
culate all  critical  orientations.  For  each  wall  edge  W  we  denote  its  two  end- 
points  as  Wi  and  W,^ ,  so  that  the  wall  region  lies  to  the  left  of  the  oriented 
segment  WiWr.  Similarly,  for  each  side  5  of  B,  we  denote  its  endpoints  by 
Si  and  Sr  ,  so  that  B  lies  to  the  iighl  of  the  oriented  segment  S^Sr  . 

Calculation  of  type  (i)  critical  orientations. 

This  phase  extends  and  provides  more  detail  of  the  procedure  briefly 
outlined  in  [LS3].  We  recall  some  definitions  and  lemmas  from  [LS3]  to  be 
used  in  the  description  of  this  phase. 

Defmtion  3.1  [LS3]:  (a)  The  tangent  line  T  of  a  contact  pair  O  =  {W,S)  of 
type  I  or  of  type  n  is  either  the  line  passing  through  W  if  Vl^  is  a  wall  edge  or 
the  line  passing  through  W  and  parallel  to  5  if  W  is  a  wall  comer  (note  that  in 
the  second  case  T  depends  on  the  orinetation  of  B). 

(b)  Let  Oi,  O2  be  two  contacts  pairs  for  which  there  exists  a  (unique)  place- 
ment Z  =  (X,  6)  of  B  at  which  it  makes  simultaneously  two  obstacle  contacts 
involving  O^  and  O2  respectively.  We  say  that  O2  bounds  Oi  at  the  orienta- 
tion e  if  fi*  =  conv  (5iU5'2)  always  intersects  W2  as  we  translate  B  from  Z 
along  the  tangent  T^  in  the  direction  of  the  intersection  point  Z12  of  Ti  and 
T2,  until  the  last  placement  at  which  Si  still  touches  W^. 

(c)  Let  Oi  =  (Wi,5i)  be  a  contact  pair.  If  Wi  is  a  wall  edge  (and  ^i  is  a 
comer  of  B),  the  two  contacts  of  ^i  against  the  two  endpoints  of  Wi  will  be 
called  endpoints  of  (the  edge  induced  by)  Oi.  If  Wi  is  a  wall  comer,  we 
define  the  endpoints  of  (the  edge  induced  by)  Oi  in  a  symmetric  manner  as 
contacts  of  the  endpoints  of  5i  against  W^. 

(d)  We  will  say  that  at  orientation  9  O2  bounds  O^  towards  an  endpoint  E  of 
the  edge  induced  by  Oi,  if  the  translation  of  B  from  Z  as  in  (b)  above  ter- 
minates at  the  contact  E. 

Proposition  3.1  [LS3]:  If  at  some  placement  Z  =  (X,  0),  B  makes  simultane- 
ously two  obstacles  contacts  involving  the  contact  pairs  Oi,  O2,  then  either 
Oi  bounds  O2  at  6,  or  O2  bounds  Oi  at  8,  except  in  the  degenerate  case 
where  the  tangents  Ti,  T2  are  coincident  or  parallel. 
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Let  Oi  be  any  contact  pair  and  let  E  be  one  of  its  endpoints.  Consider  all 
contact  pairs  that  bound  O^  (at  any  orientation  0)  towards  E.  For  each  such 
contact  pair  O2  we  define  a  bounding  Junction  FojOjC^)*  ^^^^  *^^  domain  of 
orientations  0  of  fi  in  which  O2  bounds  Oi  towards  E,  to  be  the  distance  of 
the  endpoint  of  W^i  at  which  the  contact  E  is  made  from  the  contact  point  of 
Si  and  Wi,  if  Oi  is  a  contact  pair  of  type  I,  or  tne  distance  of  the  endpoint  of 
5i  involved  in  E  from  the  contact  point  of  Si  and  Wi,  if  C>i  is  a  contact  pair 
of  type  n. 

We  partition  the  collection  of  bounding  functions  {Fo^oi)  of  Oi  into  two 
classes  Al{Oi),  A.r{Oi)  so  that  for  all  functions  FoyOi  in  ^l,  Oi  bounds  Oi 
towards  the  endpoint  Wi  of  Wi  (or  5/.  of  5i)  whereas  for  all  functions  in  A/j , 
Oi  bounds  Oj  towards  the  other  endpoint  of  Wi  (or  of  5 1).  Note  that  there 
may  exist  a  contact  pair  0/  for  which  a  bounding  function  F^jo,  appears  in 
both  collections  A/,,  A^j,  but  then  these  two  functions  will  have  openly- 
disjoint  domains  of  definition  [LS3].  Also  the  domain  of  definition  of  each 
Fq^Oi  consists  of  at  most  five  connected  intervals  [LS3].  K  the  domain  of 
such  a  Fo^Oi  ^^  ^'^^  connected,  we  regard  this  function  as  several  (^5)  distinct 
partially  defined  functions  having  connected  (and  pairwise  openly-disjoint) 
domains  (cf.  [LS3]). 

We  are  now  in  position  to  describe  the  collection  of  procedures  which 
compute  all  the  type  (i)  critical  orientations. 

Step  1:  Find  all  bounding  functions. 

For  every  two  contact  pairs  0/  =  (Wi,Si),  Oj  =  {Wj,Sj),  each  being 
either  of  type  I  or  of  type  II,  find  the  range  of  orientations  (which  can  be 
empty)  in  which  Oi  bounds  Oj  towards  a  specific  endpoint  E  of  Oj.  Split  the 
resulting  bounding  function  Fq  o^  into  (at  most  five)  "subf unctions",  and  add 
these  functions,  together  with  additional  information  concerning  their  param- 
eters and  domains  of  definition,  to  the  appropriate  collection  Ai{Oj)  or 
Ak(Oj).  (As  mentioned  above,  for  certain  pairs  Oi,  Oj  of  contact  pairs,  func- 
tions Fq  oi  may  be  added  to  both  Xi{Oj),  \ji{Oj).) 

The  results  of  [LS3]  reviewed  above  clearly  imply  that  the  storage 
required  up  to  this  point  does  not  exceed  0(k^n^).  Each  bounding  function 
(for  two  fixed  contact  pairs)  can  be  processed  in  0(1)  time,  and  thus  Step  1 
can  be  performed  in  time  O(k^n^). 

Our  next  step  is  to  find,  for  each  of  the  0{kn)  contacts  pairs  O,  the 
lower  envelope 

^L,o(Q)  =  min  {  Foo,(e)  :  Foo,  ^  A^(0)}  , 
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and  the  corresponding  envelope  '^^r  .^  for  the  collection  Ar(0). 

It  was  shown  in  [LS3]  that  etch  lype  (i)  critical  orientation  9  is  one  of 
the  three  following  kinds: 

(1)  6  is  the  orientation  of  an  intersection  point  of  two  functions  Fooi  and 
Foo  lying  along  some  lower  envelope  "^e^- 

(2)  0  is  the  orientation  of  an  intersection  point  of  some  function  Fooi  on 
^i  -o  and  a  transformation  of  a  function  Fqoj  lying  along  ^r  .q  (P'ooj  is 
shifted  and  reflected  so  that  both  fimctions  measure  distance  from  the 
same  endpoint  of  O). 

(3)  There  exist  three  contact  pairs  O^,  Oi,  O3,  such  that  0\  bounds  O2 
(towards  an  endpoint  Ei),  O2  bounds  O3  (towards  an  endpoint  £3)  and 
O3  bounds  0\  (towards  an  endpoint  £1)  at  6,  and  the  envelopes  '^e\,0\  . 
^£2;02.  ^E^fiy  ^TC  attained  at  6  by  the  functions  F03O1,  ^0i02'  -^^0203 
respectively. 

Step  2:  Calculate  lower  envelopes. 

The  following  algorithm  for  calculating  the  lower  envelopes  is  adapted 
from  [At].  Since  both  kinds  of  envelopes  are  calculated  in  much  the  same 
way,  we  will  describe  below  only  calculation  of  the  lower  envelopes  "^i  .p  for 
the  collections  ki{0). 

(1)  Fix  a  contact  pair  O,  and  partition  ^liO)  into  two  disjoint  subsets  A^,, 
A^.   of  roughly  equal  size. 

(2)  Compute  recursively  the  two  lower  envelopes 

>I^'(0)  =  min{Foo,(e):Foo,  ^  A^} 
>I^"(0)  =  min{Foo,(e):Foo,  €  A^'}  . 

Each  of  these  recursive  calculations  will  have  produced  a  sequence  of  angular 
intervals  in  each  of  which  the  corresponding  partial  envelope  is  attained  by  a 
single  bounding  function;  we  refer  to  the  endpoints  of  these  intervals  as 
breakpoints  along  the  corresponding  envelope. 

(3)  Merge  these  two  sequences  of  intervals  to  obtain  a  refined  sequence  T  of 
angular  intervals.  For  each  refined  interval  /  €  F  there  exist  unique  contact 
pairs  O',  O"  with  Fqo-  ^  ^l,  Fqo"  ^  a[' ,  such  that  ^'(0)  =  Foo'i^), 
^"(0)  =  Foo"(.Q)  for  each  0€/.  By  the  analysis  of  [LS3]  the  two  functions 
Foo'  and  Fqo"  intersect  in  at  most  4  points  (some  of  which  may  not  belong 
to  /),  which  can  be  calculated  as  the  roots  of  some  quartic  polynomial.  Each 
of  these  intersections  which  lies  in  /  is  clearly  a  breakpoint  of  "4^  =  "^ifiy  and 
every  breakpoint  of  "^  is  either  of  this  kind,  or  is  a  breakpoint  of  "*!"  or  of 
'^".   Thus  we   can   calculate   '^   from   the   merged   sequence   F   in   time 
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proportional  to  the  number  of  intervals  in  F,  which,  by  [LS3],  is  0{K6{kn)). 
Thus  the  total  time  required  by  Step  2,  when  icpeated  for  each  contact  pair  O 
and  extended  to  calculate  also  lower  envelopes  of  the  form  "^r^,  is 
0{knX6{kn)\ogkn). 

Remark:  At  the  bottom  level  of  the  recursion,  when  the  collection  of  bound- 
ing functions  consists  of  just  a  single  function,  the  sequence  of  breakpoints 
that  has  to  be  generated  simply  consists  of  the  two  endpoints  of  the  domain 
of  definition  of  that  function.  Thus,  some  of  the  breakpoints  in  the  final 
envelope  ^  will  represent  endpoints  of  such  domains. 

The  output  of  this  step  consists  of  0(kn)  sorted  Usts  of  critical  orienta- 
tions, where  each  list  represents  the  breakpoints  along  a  single  lower 
envelope  ^^.q  or  ^r.^o-  We  denote  by  ^l{0)  (rssp.  ^r{0))  the  list 
representing  "^l-o  (resp.  "^r-^o)  ^ot  a  contact  O.  For  each  such  list  ^e(0), 
the  procedure  also  produces  an  associated  list  Me{0)  of  the  bounding  func- 
tions attaining  the  corresponding  envelope  on  the  sequence  of  intervals  delim- 
ited by  the  orientations  in  ^^(O). 

By  now  we  have  found  (a  superset  of)  all  critical  orientations  9  of  type 
(i)  in  which  B  can  make  simultaneously  (at  a  free  placement)  obstacle  con- 
tacts involving  three  distinct  contact  pairs  Oi,  O2,  O3,  such  that  two  of  these 
pairs,  say  O2  and  O^,  bound  the  third  one  Oi  towards  the  same  endpoint  of 
Oi,  i.e.  such  that  the  bounding  functions  FqiOj  ^^^  ^0\02  t>elong  to  the 
same  collection  Ai,(C>i)  or  kR{0{)  and  meet  one  another  at  9  along  the 
corresponding  lower  envelope.  (Note  that  not  all  orientations  corresponding 
to  envelope  breakpoints  necessarily  yield  free  critical  contacts  of  B.  How- 
ever, subsequent  steps  of  the  algorithm  will  be  able  to  detect  such  spurious 
breakpoints  and  discard  them;  see  below  for  more  details.) 

Step  3:  Calculate  the  second  kind  of  type  (!)  critical  orientations. 

These  are  orientations  9  at  which  B  makes  simultaneously,  at  some  free 
placement,  obstacle  contacts  involving  three  distinct  contact  pairs  Oi,  C>2,  O3, 
such  that  two  of  them,  say  O2  and  O-^,  bound  O^,  but  with  F^jOj  belonging 
to  A^(Oi)  and  F01O3  belonging  to  A/j(C>i)  at  9.  Calculation  of  such  orienta- 
tions is  easy  and  proceeds  as  follows. 

(1)  For  each  contact  pair  O,  merge  the  two  hsts  O^.  {O)  and  <l>^  {O)  of  break- 
points to  obtain  a  refined  sorted  list  of  angular  intervals  delimited  by  these 
breakpoints. 

(2)  For  each  such  refined  interval  /,  there  exist  unique  functions 
^00 1  ^  ^l{0),  F002  ^  A^(0)  which  attain  ^i.q  and  ^r.q  respectively  over 
/.  Calculate  the  (at  most  4)  intersections  between  these  two  functions  and  add 
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those  intersections  that  he  in  /  to  the  output  Hst  of  new  critical  orientations. 
(Note  that  in  this  calculation  t'ooi  ^^  ^°  ^^  shifted  and  reflected  into  the 
coordinate  frame  in  which  Fqoi  i»  defined,  so  that  both  functions  will  meas- 
ure the  distance  from  the  same  endpoint  of  O.) 
Qearly  Step  3  nms  in  0(kn\(,{kn))  time. 

Step  4:  Calculate  the  third  kind  of  type  (i)  critical  orientations. 

Finally  we  calculate  the  third  and  most  complex  kind  of  critical  orienta- 
tions. At  each  such  orientation  8  B  can  make  simultaneously  a  free  triple 
contact  involving  three  distina  contact  pairs  01,02,0^,  such  that 
F0102  ^  A£j((9i),  F0202  ^  ^£2(02),  Fo^oi  ^  AfjCOs),  where  E,  ^  {L,R} 
for  1  =  1,2,3  (and  such  that  ail  thefie  functions  lie  on  the  corresponding  lower 
envelopes).  In  other  words,  no  two  of  these  contact  pairs  simultaneously 
bound  the  third  one. 

To  find  these  orientations,  we  first  merge  all  Hsts  ^l(0),  0^(0),  for  all 
contact  pairs  O,  to  obtain  a  single  sorted  list  4>  consisting  of  0{knXs{^n)) 
refined  noncritical  intervals  delimited  by  all  these  breakpoints.  Each  interval 
/  in  O  has  the  property  that  each  envelope  '^e;o  is  attained  over  /  by  a  single 
bounding  function  Fqo'  • 

Step  4  processes  the  intervals  in  0  in  order,  maintaining  a  priority  queue 
Q  of  potential  critical  orientations  of  the  kind  we  seek,  and  a  map  ^*  which 
maps  each  pair  {0,E)  of  a  contact  pair  O  and  endpoint  E  to  the  contact  pair 
O'  for  which  "^e.o—P^oo'  over  the  currently  processed  interval  of  <i>.  This 
procedore  runs  as  follows. 

(1)  Starting  at  the  initial  interval  /q,  we  first  initialize  the  value  of  "^^  for 
/o,  and  then  construct  a  list  A  of  potential  triple  contacts  of  the  third  kind 
induced  by  the  various  lower  envelopes  over  /q,  as  follows.  For  each  contact 
pair  Oi,  and  each  E^  ^  {L,R},  find  the  (unique)  contact  pair  O2  such  that 
'^Ei;Oi=P'oi02  over  /o  (using  the  map  ^*).  For  each  £2  ^  {^.^}  find  the 
(unique)  contact  pair  C>3  such  that  '^e2;02'^Fo203  over  Iq.  Finally  for  each 
£3  €  {L,R]  for  which  ^Ey,02,-Fo-iOi  over  /q,  add  the  item 
{Oi,02,0-i,Ei,E2,E2,)  to  A.   Qearly,  A  consists  of  at  most  0{kn)  items. 

(2)  Next  initialize  the  priority  queue  Q  as  follows.  For  each  item 
{0\,02,02,,Ei,E2,E2,)  in  A,  calculate  the  (at  most  4)  orientations  at  which  a 
triple  contact  involving  the  pairs  0^,02,02,  can  be  made  such  that  at  this 
contact  O2  bounds  0^  towards  £1,  Oj,  bounds  O2  towards  E2,  and  Oi 
bounds  O3  towards  £3.  Orientations  of  this  kind  that  also  belong  to  /q  are 
immediately  added  to  the  output  list  of  critical  orientations,  whereas  all  the 
others  are  added  to  Q,  each  tagged  by  the  corresponding  item  (Oi,  .  .  .  ,£3). 
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(3)  Next  we  process  the  remaining  intervals  in  O  in  order.  Let  /  be  the 
interval  presently  being  processed,  let  /'  be  the  preceding  interval,  and  let  6 
be  the  common  endpoint  of  /'  and  /.  By  our  assumptions  at  most  one  triple 
contact  can  occur  at  6.  Suppose  this  contact  involves  three  individual  contact 
pairs  Oi, 02,03  such  that,  say,  O2  and  O3  both  bound  Oi  towards  its  L- 
endpoint.  Then  the  only  envelope  for  which  the  bounding  function  attaining 
it  changes  from  /'  to  /  is  "^l-.Oi  •  Thus  if  we  were  to  calculate  (in  the  manner 
described  above)  a  list  A  of  potential  triple  contacts  from  the  envelopes  over 
/,  we  would  have  obtained  the  same  list  that  would  be  obtained  from  the 
envelopes  over  /',  with  the  exception  of  contacts  involving  Oi  (more  pre- 
cisely, contacts  corresponding  to  items  of  the  form  (Oi,0' ,0" ,L,E' ,E")). 
But  there  are  clearly  only  0(1)  items  of  this  form,  and  they  can  all  be  com- 
puted in  0(1)  time  using  the  map  '^*  (for  the  interval  I)  as  in  (1)  above. 

We  therefore  proceed  as  follows. 

(3.i)  Update  the  value  of  ^*{Oi,L),  obtain  the  0(1)  new  potential  triple 
contacts  arising  from  this  change,  and  insert  the  corresponding  critical  orien- 
tations into  Q.  (Optionally,  we  may  also  choose  to  delete  the  0(1)  entries 
from  Q  that  were  obtained  using  the  old  value  of  ^*{Oi,L).) 

(3.ii)  Remove  from  Q  all  entries  whose  orientations  He  in  /.  Let 
(e,(0,0',0 ",£,£',£:"))  be  an  entry  removed  from  Q.  We  check  whether 
this  contact  is  still  valid,  in  the  sense  that  the  present  value  of  ■4''*  satisfies 
^*(0,£:)  =  0',  ^*{0',E')=0",  and  >I'*(0",£")  =  0.  If  so,  add  this  entry 
to  the  output  list. 

This  concludes  the  description  of  Step  4.  Its  time  complexity  is  clearly 
0{kn\^{kn)  log  kn).  Therefore  all  critical  orientations  of  type  (i)  can  be  cal- 
culated in  time  0(knX(,{kn)  log  kn). 

Step  5:   Calculation  of  type  (ii)  critical  orientations: 

These  orientations  are  calculated  in  a  similar  (and  simpler)  manner  to  the 
calculation  of  type  (i)  orientations.  Recall  that  a  type  (ii)  orientation  6  is  one 
in  which  a  vertically  maximal  corner  u  of  Kq  induced  by  a  type  III  contact 
pair  O  =  {W,S)  lies  immediately  below  another  (convex  or  non-convex) 
comer  u*  of  that  set.  If  u*  is  convex,  then  «*  is  induced  by  another  contact 
pair  (Wi,Si)  of  a  wall  comer  W^  and  a  comer  S^  of  B.  In  this  case  6  can  be 
found  by  intersecting  the  circle  at  radius  \SSi\  about  W  with  the  vertical  ray 
descending  from  W^.  There  are  at  most  two  such  intersections,  and  for  each 
of  these  intersection  points  Z  we  obtain  a  corresponding  critical  orientation  6 , 
at  which  the  segment  55 1  on  fi  becomes  parallel  to  WZ.  It  is  therefore  clear 
that  the  total  number  of  such  critical  orientations  is  at  most  0{k^n^),  and 
that,  by  repeating  the  above  procedure  for  every  two  type  DI  contact  pairs 
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{W,S),  (Wi,5i),  we  obtain  in  time  0{k'^n'^)  a  superset  of  the  set  containing 
all  these  critical  orientations. 

K  the  comer  w*  is  non-convex,  ihen  we  have  to  apply  the  techniques 
described  in  Subsection  3.3  to  calculate  the  corresponding  critical  orienta- 
tions.  Specifically, 

(1)  For  each  type  HI  contact  pair  0=(W,S)  and  for  each  contact  pair 
0,  =  (w,,Si)  of  type  I  or  IT,  find  (in  constant  time)  the  domain  of  defini- 
tion and  other  parameters  of  the  function  Gooii^)-  (Recall  that  this 
function  measures,  for  each  orientation  0,  the  vertical  distance  between 
the  comer  u  of  ATe  induced  by  O  and  the  edge  e  of  that  set  induced  by  O,, 
whenever  this  distance  is  well-defined  and  u  lies  below  e.) 

(2)  Calculate  the  lower  envelope  '¥o  of  tiie  functions  Gooi,  over  all  contact 
pairs  Oi ,  using  the  procedure  as  applied  in  step  2  to  the  functions  Fooi  ■ 
This  produces  a  list  Tq  of  critical  orientations,  each  representing  a  break- 
point along  ^o-  By  the  results  in  Subsection  3.2,  Tq  consists  of  at  most 
0(\e{kn))  orientations. 

(3)  Finally,  we  can  merge  Tq  with  the  previously  calculated  Ust  of  potentially 
critical  orientations  9  at  which  u  ties  below  another  convex  comer  of  ATe- 
We  can  determine  whether  any  such  latter  orientation  6  is  indeed  critical 
by  testing  whether  it  represents  a  point  actually  lying  on  the  envelope 
■^Po,  and  discard  6  if  it  is  found  not  to  be  critical. 

Qearly  step  5  of  our  algorithm  requires  0(kn\6{kn)  log  kn)  time. 

Step  6:  Calculation  of  type  (ill)  critical  orientations 

Recall  that  this  type  of  critical  orientations  arises  when  two  adjacent 
edges  of  Kq  become  collinear,  or  when  an  edge  of  A^e  becomes  horizontal. 
As  noted  in  Subsection  2.2,  there  are  at  most  0{kn)  orientations  of  the  first 
kind  that  are  not  also  type  (i)  critical  orientations,  and,  assuming  general 
position  of  the  walls  (so  that  no  wall  edge  is  horizontal),  there  are  at  most 
0{k)  orientations  of  the  second  kind.  All  these  orientations  are  trivial  to  cal- 
culate in  time  0{kn)  (note  that  as  before  this  calculation  produces  a  superset 
of  the  required  critical  orientations,  because  not  all  these  orientations  neces- 
sarily arise  in  semi-free  critical  contacts  of  B). 

3.4.   Constructing  the  edge-graph  EG 

By  now  we  have  calculated  a  set  T  of  size  OiknX^ikn))  that  contains  all 
the  critical  orientations  of  all  three  types.  Suppose  that  T  is  sorted  in  ascend- 
ing order,  and  assume  without  loss  of  generality  that  0  =  0  is  not  an  orienta- 
tion in  T. 
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As  in  Section  2,  we  represent  each  node  i  of  the  graph  EG  by  a  pair 
(u,L^),  where  L„  =  (61,82)  is  the  angular  life  span  of  4,  and  where  u  is  the 
(discrete  labeling  of  the)  comer  of  VGq,  for  0  C;  L„ ,  that  lies  on  the  edge  of 
FP  represented  by  ^. 

The  calculation  of  EG  will  be  accomplished  in  a  manner  similar  to  that  of 
[LSI].  That  is,  we  process  critical  orientations  in  increasing  order,  maintain- 
ing the  "cross-section"  graph  VGq  and  use  it  to  update  EG  at  each  critical 
orientation.  At  each  such  orientation  6  we  determine  those  nodes  of  EG 
whose  life-span  terminates  at  6  (these  are  nodes  whose  corresponding  comers 
have  to  be  deleted  from  VGe),  and  the  new  nodes  whose  life  span  starts  at  9. 
Nodes  of  the  first  kind  will  already  have  been  stored  in  EG,  and  we  update 
their  life  span  by  adding  6  as  its  terminal  orientation.  Nodes  of  the  second 
kind  are  added  to  EG,  with  6  as  the  initial  orientation  of  their  life  span,  and 
with  the  corresponding  terminal  orientation  being  presently  left  undefined. 
EG  is  also  augmented  by  edges  connecting  the  nodes  just  inserted  with  nodes 
already  present  in  EG,  as  described  in  Subsection  2.2  (note  that  these  edges 
correspond  to  new  edges  in  VGe' ,  for  0'  slightly  larger  than  9). 

To  initialize  EG,  we  first  calculate  the  graph  VGe=o-  For  each  comer  u 
of  this  graph  we  create  a  corresponding  initial  node  in  EG  whose  label  con- 
tains u,  but  with  both  initial  and  terminal  life-span  orientations  left  unde- 
fined. 

Our  algorithm  has  thus  the  following  structure: 

(1)  Merge  all  the  various  lists  of  critical  orientations  into  a  single  sorted  hst 
T  consisting  of  0  <  Bi  <  62  <  •  •  •  <  9^  <  27t  (where 
w  =  0(kn\(,(kn))).  (Without  loss  of  generality  we  assume  that  0  is  not 
a  critical  orientation.) 

(2)  Compute  VGq  for  9  =  0.  Copy  VGe  into  EG,  expanding  the  label  of  each 
corner  and  edge  of  VG  e  by  two  (presently  undefined)  entries  for  the  two 
delimiting  orientations  of  the  life-span  of  the  resulting  node  or  edge  of 
EG. 

(3)  For  each  critical  orientation  9/,  i  =  1,  .  .  .  ,w,  do: 

(a)  Determine  the  changes  that  VG^  undergoes  as  9  crosses  9;  from 
smaller  to  larger  orientations. 

(b)  Let  V 1 ,  .  .  .  , v/  be  the  comers  and  ei,  .  .  .  ,e„  he  the  edges  of  VG e 
that  disappear  as  9  crosses  9/.  Remove  them  from  VGq,  and  assign 
9/  as  the  terminal  life-span  orientation  of  the  corresponding  nodes 
and  edges  in  EG. 

(c)  Let  v*  ,  .  .  .  ,v*  be  the  comers  newly  appearing  in  VGq  after  9 
crosses  9,,  and  let  e*  ,  .  .  .  ,e*  be  the  new  edges  in  that  graph.  Add 


-28- 

these  elements  to  VG  9 ,  and  add  the  corresponding  nodes  and  edges 
to  EG  with  6/  as  their  iridal  life-span  orientation  (and  with  their  ter- 
minal life  span  orientation  presently  left  undefined). 

(4)  {At  the  end  of  the  above  main  loop  the  nodes  and  edges  of  EG  whose 
terminal  life-span  orientation  is  still  undefined  are  in  1-1  correspondence 
with  the  initial  nodes  and  edges  of  EG,  as  created  in  step  (2)  above  from 
the  initial  graph  VGe=o-} 

This  final  step  in  the  construction  of  EG  simply  identifies  these  two  types 
of  nodes.  Specifically,  for  each  node  (w,(9i,fl))  of  EG  whose  terminal 
life-span  orientation  is  still  undefined,  find  the  corresponding  node 
(M,(n,02))  of  EG  whose  initial  life-span  orientation  is  undefined  and 
merge  them  into  the  single  node  (",(01,62)).  (Note  that  there  are  only 
0{kn)  such  pairs  of  nodes.) 

Remark:  In  step  3(a)  we  can  first  check  whether  the  criticality  that  is  sup- 
posed to  happen  at  6  actually  arises  in  VG  e ;  otherwise  6  is  a  spurious  critical 
orientation  which  we  can  simply  discard.  (For  example,  if  0  is  an  orientation 
at  which  B  makes  three  obstacle  contacts  simultaneously,  then,  if  this  is  really 
an  obstacle  contact  then  for  6'  near  6  the  boundary  of  ATe'  must  contain  edges 
induced  by  all  these  three  contacts,  and  we  can  easily  check  whether  this  is 
the  case  from  the  current  value  of  VG  e .) 

We  can  now  state  our  main  result. 

Theorem  3.2:  The  graph  EG  has  at  most  0{kn\e(kn))  nodes  and  edges.  The 
algorithm  just  described  computes  EG  in  0{kn\^{kn)  log  kn)  time. 

Proof:  Let  0o  be  a  critical  orientation  in  T,  and  let  (Oq  — e,0o  +  €)  be  a  suffi- 
ciently small  interval  in  which  only  0o  is  critical.  As  in  Section  2,  we  also 
make  the  simplifying  assumption  that  only  one  critical  contact  of  B  occurs  at 
00.  (See  also  a  discussion  on  this  assumption  in  [LS3].  An  exception  are 
orientations  of  type  (iii)(b),  at  each  of  which  up  to  0{n)  different  expanded 
obstacle  edges  become  simultaneously  horizontal;  see  below  for  the  special 
treatment  of  these  orientations.)  The  changes  that  occur  in  VGq  (and  induce 
similar  changes  in  EG)  as  6  varies  from  Oq-c  to  ©o  +  e  depend  on  the  type  of 
criticality  that  arises  at  0o : 

(a)  Suppose  first  that  0o  is  a  critical  orientation  of  type  (i)  at  which  three 
edges  of  Kq  become  concurrent  at  some  non-convex  comer  of  Kq.  Then 
VGq  can  change  in  one  of  the  following  ways:  either 

(i)  one  corner  of  this  graph  disappears  and  is  replaced  by  two  new 
corners  connected  by  a  new  edge  (see  Fig.  3.2(a));  or 
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(ii)  two  adjacent  comers  of  VG  e  merge  and  fonn  a  new  comer,  and  the 
edge  connecting  them  disappears  (see  Fig;  3.2(a));  or 
(iii)  a  small  connected  component  of  VG  e  consisting  of  just  three  comers 
shrinks  to  a  single  isolated  corner  and  then  disappears;  (see  Fig.  3.2(b)) 
or,  finally, 

(iv)  a  new  connected  component  of  VGq  newly  appears,  first  as  an  iso- 
lated comer  and  then  expands  to  a  triangle  (see  Fig.  3.2(b)). 

Since,  by  the  arguments  in  Subsection  3.1,  the  degree  of  each  vertex  of 
VG  e  at  any  noncritical  6  is  at  most  three,  it  follows  that  in  each  of  these  cases 
step  (3)  of  the  algorithm  newly  generates  or  modifies  only  a  small  fixed 
number  of  nodes  and  edges  of  EG  as  a  result  of  these  changes  in  VG^  at 
e  =  Bo. 

(b)  Next  suppose  that  8o  is  a  critical  orientation  of  type  (i)  at  which  a  convex 
comer  of  Kq  comes  to  lie  on  an  edge  of  Kq.  In  this  case  VG^  can  change  in 
one  of  the  following  ways:  either 

(i)  one  convex  comer  of  that  graph  disappears  and  is  replaced  by  two 
new  nonconvex  corners;  these  comers  lie  in  the  interior  of  a  former  edge 
of  VGq,  and  now  split  that  edge  into  two  new  edges  of  that  graph,  each 
dehmited  by  one  of  these  new  corners  (see  Fig.  3.3(a));  or,  conversely, 

(ii)  two  nonconvex  corners  u,  v  of  VGe  merge  into  a  new  convex  corner 
w;  furthermore,  w  becomes  adjacent  to  one  neighbor  u'  of  u  and  to  one 
neighbor  v'  of  v,  and  the  two  remaining  neighbors  s,t  (of  u  and  v  resp.) 
become  adjacent  to  one  another  (see  Fig.  3.3(a)). 

Thus  in  each  of  these  cases  EG  undergoes  only  0(1)  changes  at  Qq. 

Since  there  are  at  most  0{kn\(,{kn))  critical  orientations  of  type  (i),  the 
total  number  of  new  nodes  and  edges  added  to  EG  at  these  orientations  is 
also  bounded  by  0{kn\(,{kn)). 

(c)  Next  suppose  that  0o  is  a  critical  orientation  of  type  (ii),  in  which  a  verti- 
cally maximal  corner  u  of  Kq  comes  to  lie  directly  below  another  comer  v  of 
ATe  (i.e.  the  associated  auxiliary  comer  u*  coincides  with  v).  In  this  case  u* 
changes  its  location  from  one  edge  (or  portion  of  an  edge)  vw^  of  ATe  to  an 
adjacent  one  vh'2.  Consequently  the  two  former  edges  vu* ,  u*wi  of  VGq 
have  to  be  merged  into  a  single  edge  vwi,  whereas  the  edge  vh'2  is  split  into 
two  new  edges  vw*,  u*w2  (See  Fig.  3.3(b);  note  that  the  new  edge  vu*  is  a 
portion  of  an  expanded  obstacle  edge  induced  by  a  contact  pair  different  from 
that  inducing  the  previous  edge  vu*).  Again  only  0(l)  changes  have  to  be 
made  at  6o  in  the  graphs  VGq  and  EG. 

(d)  Finally  suppose  6o  is  a  critical  orientation  of  type  (iii).   Suppose  first  that 
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8o  is  of  subtype  (iii)(a),  i.e.  a  side  5  of  B  at  orientation  80  becomes  parallel 
to  a  wall  edge  W.  Let  Wi ,  W^  L  e  the  endpoints  of  W  and  let  S^ ,  Sr  be  the 
endpoints  of  5.  Without  loss  of  generality  suppose  that  for  6  slightly  less 
than  60,  ATe  has  two  adjacent  edges  induced  respectively  by  the  contact  pairs 
(WtyS),  (W,Sl),  and  that  these  edges  become  collinear  at  Gq,  and,  for  9 
sHghtly  larger  than  0o  are  replaced  by  two  other  edges  induced  respectively 
bythepairs(W,5^),(M^^,5). 

Consequently,  for  arbitrarily  small  e>0,  all  the  vertices  and  edges  of 
FGe-t  that  include  (Wi,S)  or  (W,Si)  in  their  label  terminate  at  9,  and 
corresponding  new  vertices  and  edges  that  include  {W,Sj{),  {Wj(,S)  in  their 
label  show  up  in  VGe+t-  The  rest  of  the  vertices  and  edges  of  both  VG 
graphs  remain  unchanged.  Since  there  are  at  most  0(kn)  edges  and  vertices 
including  these  pairs  in  their  label  in  VGe-«  (resp.  VGe+t),  and  the  degree  of 
each  of  these  vertices  is  at  most  three,  the  number  of  edges  and  vertices 
removed  from  (resp.  added  to)  VGq  and  EG  at  each  critical  orientation  of 
type  (iii)(a)  is  at  most  0{kn).  Since  there  are  at  most  0{kn)  critical  orienta- 
tions of  this  type,  the  total  number  of  modifications  of  VGq  and  of  EG  at 
these  orientations  is  at  most  0{k^n^). 

Next  suppose  that  60  is  of  subtype  (iii)(b),  i.e.  a  side  S  of  B  becomes 
horizontal  (with  B  lying  above  S)  at  9o .  Let  5^ ,  Sr  be  the  endpoints  of  S. 
Then  without  loss  of  generality  we  can  assume  that  for  9  slightly  less  (resp. 
larger)  than  Gq,  5/.  is  lower  (resp.  higher)  than  Sr.  Then  for  9  slightly  less 
than  9o,  Kq  contains  at  most  0(n)  vertically  maximal  comers  induced  by  con- 
tact pairs  of  the  form  (W,Si),  which  cease  to  be  vertically  maximal  as  Gq  is 
crossed,  and  have  to  be  replaced  by  another  set  of  at  most  0(n)  adjacent 
comers  induced  by  corresponding  pairs  of  the  form  (W,Sr)  which  newly 
become  vertically  maximal.  Thus,  arguing  as  above,  at  most  0(n)  changes 
are  required  in  the  graphs  VG^  and  EG  at  each  of  the  0(k)  critical  orienta- 
tions of  type  (iii)(b),  thus  the  total  number  of  graph  modifications  at  these 
orientations  is  only  0{kn). 

This  completes  the  proof  of  the  first  part  of  Theorem  3.2. 

Analysis  of  the  time  complexity  of  the  algorithm:  We  have  already  shown 
that  the  preliminary  stages  of  our  procedure,  which  calculate  all  critical  orien- 
tations, require  0{kn\^{kn)  log  kn)  time.  As  to  the  complexity  of  the  main 
part  of  the  algorithm  described  in  this  Subsection  we  argue  as  follows. 

In  Step  (1)  we  merge  0{kn)  sequences  of  critical  orientations,  each  hav- 
ing length  at  most  0{\(,{kn));  this  is  easily  accomplished  in  time 
0{kn\(,{kn)\ogkn). 

Construction  of  VGe  =  o  in  Step  (2)  can  be  accomplished  e.g.  in  time 
0{kn  \o^kn)  using  the  algorithm  described  by  [KS]. 
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As  to  Step  (3),  it  processes  0{kn\e{f^n))  critical  orientations.  At  each 
orientation  6,  of  type  (i)  or  (ii)  it  has  to  perfonrji  oi;ly  0(1)  updating  opera- 
tions on  the  graphs  VG^  and  EG,  each  of  wliich  gan  be  accomplished  in  time 
0{\ogkn).  ■■'■ 

If  6/  is  a  critical  orientation  of  type  (iii)(a),  let  ej,  ei  be  the  two  adjacent 
expanded  obstacle  edges  that  become  collinear  at  9/  and  are  then  replaced  by 
two  other  edges,  e^,  e^  (see  the  preceding  paragraph  (d)  in  this  proof). 
Then  it  is  easily  checked  that  the  boundary  of  FP^  varies  continuously  with  0 
(in  the  Hausdorff  topology  of  sets),  and  also  all  the  vertices  of  VGe  vary  con- 
tinuously with  e  near  9/.  It  follows  that  the  combinatorial  changes  that  VGq 
undergoes  at  9;  are  just  changes  of  the  labels  of  all  vertices  of  VG^  that  lie  on 
e-i  U  62  before  crossing  9/  (and  become  points  on  ^3  U  e^  after  crossing  9/), 
and  of  their  incident  edges.  Since  VG^  contains  at  most  0(kn)  such  vertices 
and  edges,  and  since  they  can  all  be  accessed  in  total  0{kn  log  kn)  time  using 
our  representation  of  this  graph,  we  can  update  VGq  (and  EG)  during  pro- 
cessing of  the  0{kn)  orientations  of  type  (iii)(a)  in  overall  time 
0(k^Ti^  log  kn). 

Finally,  if  9,  is  a  critical  orientation  of  type  (iii)(b),  let  5  be  the  side  of  B 
that  becomes  horizontal  at  9/,  and  let  S^,  Sj^  be  its  endpoints.  As  in  the 
above  discussion,  to  update  VGq  at  9^,  one  has  to  remove  from  this  graph  the 
auxiliary  comers  corresponding  to  all  previous  vertically  maximal  comers 
induced  by  contact  pairs  of  the  form  {W,Si),  and  also  remove  all  their 
incident  edges,  and  then  add  new  auxiliary  corners  corresponding  to  new 
vertically  maximal  comers  induced  by  pairs  of  the  form  (W,5/j),  and  also  add 
their  incident  edges.  Since  these  new  auxiliary  comers  are  not  yet  present  in 
VGe,  we  need  to  create  them  and  locate  them  on  the  boundary  of  Kq  in  order 
to  determine  their  adjacent  vertices  along  this  boundary.  This  can  be  accom- 
plished by  a  standard  and  straightforward  line-sweeping  technique,  running  in 
time  0(kn  log  kn).  Since  the  number  of  type  (iii)(b)  orientations  is  only 
0{k),  processing  of  these  orientations  at  Step  (3)  of  the  algorithm  requires 
only  0{k^n  log  kn)  time. 

In  the  final  Step  (4),  all  we  need  to  do  is  to  pair  each  vertex  ^  = 
(M,(n,9/))  of  EG  whose  initial  life-span  orientation  is  undefined,  with  a 
corresponding  vertex  ^'  =  (u,  (9i,n))  having  the  same  label  u  and  an  unde- 
fined terminal  life-span  orientation.  For  each  such  pair  we  form  a  new  ver- 
tex i*  =  (u,  (9i,92))  of  EG  and  merge  the  two  lists  of  edges  incident  to  C,  C' 
to  a  single  list  for  ^*.  Since  there  are  only  0(kn)  such  pairs  of  vertices,  this 
merging  can  be  trivially  accomplished  in  0{kn)  time.  Hence  the  time 
required  by  the  algorithm  is  at  most  0(knX(,(kn)  log  kn),  as  asserted.  This 
completes  the  proof  of  the  theorem.  □ 
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Remarks:  (1)  As  noted  earlier,  when  m<<n  the  complexity  of  VG^  for 
any  fixed  0  is  only  0{km''n)  which  can  be  significantly  smaller  than  0{kn). 
We  can  exploit  this  fact  to  rtbt'Jn  a  better  bound  on  the  complexity  of  our 
algorithm  as  follows. 

For  each  of  the  0{kn)  contact  p^iifi-  O  =  iW,S)  let  Jq  denote  the  angular 
interval  of  orientations  9  of  B  at  which  an  obstacle  contact  induced  by  O  is 
possible  (in  the  sense  of  Definition  2.1(c)).  The  0(kn)  endpoints  of  these 
intervals  partition  the  angular  .space  of  orientations  into  0{kn)  disjoint  inter- 
vals /i,/2,  .  .  .  ,/,.  It  follows  from  the  argimients  of  [GRS]  that  within  each 
of  these  intervals  Ij  the  set  of  possible  obstacle  contacts  (in  the  sense  of 
Definition  2.1(c))  remains  constant,  and  that  there  are  at  most  km+n  such 
contacts.  Hence  the  number  of  pairs  of  contact  pairs  that  can  induce  a  double 
contact  of  B  with  an  orientation  In  Ij  is  at  most    f     o^")  ^  0((km+n)^). 

Moreover,  as  we  cross  from  Ij  io  the  next  interval  Ij  +  i,  only  one  contact  pair 
can  induce  a  new  obstacle  contact,  thus  the  number  of  new  possible  pairs  of 
contacts  within  Ij  +  i  is  only  km+n,  and,  repeating  this  argument  over  all 
intervals  //  we  obtain  that  the  total  number  of  double  contacts  of  B  is  at  most 

^^"2^"]  +  (km+n)  0{kn)  =  0{{km+n)kn). 

For  each  contact  pair  O,  let  to  denote  the  number  of  bounding  functions  Fqo- 
appearing  in  the  envelopes  ^i.;o,  ^r;o-  Then  the  combinatorial  complexity 
of  ^L;o^  ^R;o  is  0(\6{to)),  SO  that  the  total  number  of  type  (i)  critical 
orientations  is  ^  ^^sih)  which,  as  can  be  easily  verified  (using  the  bound  in 

o 
[Sh]),  is  O {{km +n)\eikn)). 

Similar  arguments  apply  to  all  other  types  of  critical  orientations,  and  there- 
fore imply  that  the  complexity  of  our  algorithm  is  at  most 
0{{km+n)  \e{kn)  \og  kn). 

(2)  The  example  given  in  [LS3]  can  be  used  to  show  that  the  combina- 
torial comlexity  of  FP  can  be  Cl{k^n^)  in  the  worst  case,  so  that  our  algo- 
rithm is  close  to  being  optimal  among  algorithms  that  calculate  the  entire 
space  FP.  See  also  a  related  result  by  O'Rourke  [OR]  for  the  case  of  a  mov- 
ing line  segment. 

3.5.   Actual  motion  planning. 

Once  EG  is  available,  actual  motion-planning  between  any  two  given 
placements  of  B  can  be  easily  accomplished,  as  follows.  Let  Z  =  (X,  6), 
Z'  =  (X'  ,6')  be  two  given  free  placements  of  B.  We  first  calculate  the  nodes 
0(Z),  0(Z')  of  EG  (where  4>  is  the  map  defined  in  the  paragraph  following 
the  proof  of  Proposition  2.6).    It  is  easy  to  calculate  these  nodes  in  0{kn) 
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time  in  a  straightforward  manner  following  the  definition  of  <I>.  Next  deter- 
mine by  simple  graph  searching  whether  0(Z)  and  0(Z')  belong  to  the  same 
connected  component  of  EG  (this  can  be  acfccmplished  iii  time  O (knkeil^n))) . 
If  not,  then  by  Theorem  2.2  no  collision-free  motion  of  B  between  place- 
ments Z  and  Z'  is  possible.  •  Otherwise,  let 
IT  =  (Ci  =  <I>(Z),^2,  ■  ■  '  ,it  =  ^{Z') )  be  a  path  in  EG  connecting  0(Z)  to 
<I>(Z').  We  transform  tt  into  a  continuous  semi -free  motion  of  B  from  Z  to 
Z'  as  follows.  First  B  is  translated  from  Z  to  a  placement  on  the  edge  ^i  of 
FP,  as  in  the  definition  of  the  map  O.  Then,  for  each  /=1,  .  .  .  ,f-l  we 
choose  some  orientation  0,  in  the  life  span  of  the  edge  (^/,^/  +  i)  of  EG,  and 
move  B  from  its  current  position  on  ^;  to  a  position  on  ^/  +  i  in  two  substeps: 

(a)  First  B  is  moved  along  the  edge  ^/  until  its  orientation  becomes  equal  to 
0/.  The  exact  nature  of  this  motion  depends  on  the  type  of  comer  u  of 
ATe  that  labels  |,.  For  example,  if  w  is  a  convex  comer  of  Kq,  induced  by 
a  contact  pair  (W,5)  of  type  III,  then  the  required  motion  of  B  is  simply 
a  rotation  of  B  about  its  comer  S  which  remains  touching  the  wall  corner 
W.  If  w  is  a  non-convex  comer  of  Kq,  induced  by  two  contact  pairs 
(^i^'S'i),  C^^2.'S'2)  of  type  I  or  n,  5  has  to  glide  as  it  maintains  the  two 
contacts  of  ^i  against  Wi  and  of  S2  against  W2.  The  resulting  motion, 
known  as  a  glissette,  can  be  parametrized  by  a  quartic  polynomial  equa- 
tion (see  [SS]).  Finally,  if  «  is  an  auxihary  comer  v*  corresponding  to 
some  convex  comer  v  ot  K^,  let  (W,5)  be  the  contact  pair  inducing  the 
edge  of  ATe  containing  u,  and  let  {W ,S')  be  the  type  III  pair  inducing  v. 
Then  B  has  to  move  so  that  it  maintains  contact  of  S  against  W  while  its 
corner  S'  remains  co-vertical  with  (and  above)  the  wall  comer  V/' .  This 
motion  is  clearly  also  a  fourth  degree  glissette  (where  we  can  think  of  S' 
as  maintaining  contact  with  the  vertical  ray  emanating  upwards  from 
W). 

(b)  After  B  has  reached  orientation  0,,  we  next  move  it  to  the  FP-edge  4/  + 1 . 
By  definition,  VGe,  contains  an  edge  (w,v)  where  u  (resp.  v)  is  the  label 
of  the  corner  oi  K^  that  induces  4/  (resp.  C/  +  i)-  Thus  the  required 
motion  of  B  can  be  accomplished  by  translating  it  along  the  straight 
boundary  segment  of  K^.  connecting  u  and  v. 

Finally,  after  reaching  the  last  edge  ic,,,  we  complete  the  motion  of  B  by  first 
moving  it  along  ^,  until  it  reaches  its  final  orientation  0'  (as  in  (a)  above), 
and  then  translate  it  back  to  Z'  by  applying  the  reverse  of  the  canonical  trans- 
lation used  in  the  definition  of  the  map  O. 

Remark:  The  motion  of  B  described  above  is  semi-free  rather  than  free. 
This  is  a  consequence  of  our  choice  of  representation  of  FP,  and  resembles 
the  motion  obtained  for  a  moving  line  segment  [SiS]. 
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Fig.  1.1.  An  instance  of  the  motion  planning  problem 


Fig.  2.1.  ^e  and  FP( 
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Fig.  2.2:  FPe  and  VGe 
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Fig.  2.3(a) 


Fig.  2.3(b) 
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Fig.  2.3(c) 
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Fig.  2.3(e) 


39 


Fig.  2.4 
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Fig.  3.1(a) 


Fig.  3.1(d) 


41 


0, 

/ 

/ 

< 

1- 

S 

s 

02 

/ 

r^ 

N 

0 


Fig.  3.2(a) 


\ 


\ 


/ 


\ 
\ 


/' 


/ 


/ 


\ 


01 


\ 


03 

\0o 


\ 


\ 


r. 


/61 


\ 


/ 


/ 


m 


0- 


\ 


\0; 


Fig.  3.2(b) 


42 


Fig.  3.3(a) 


Fig.  3.3(b) 
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